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Abstract 


X 

TThe  Kolb-Griem  theory  of  spectral  line  broadening  occurring  in 
plasmas  is  discussed;  emphasis  is  placed  on  the  importance  of  the 
electric  microfield  distribution.  The  Mozer-Baranger  calculation  of 
the  electric  microfield  distribution  in  a  plasma  is  outlined;  its 
dependence  upon  the  nature  of  the  component  particle  correlations  is 
shown.  The  electrons  and  ions  must  be  treated  separately  because  of 
their  inherent  differences;  the  electrons  interact  through  a  Coulomb 
field,  whereas  the  ions  interact  with  shielded  fields.  The  Salpeter 
cluster-expansion  technique  is  considered  in  detail.  It  is  applied 
to  the  determination  of  pair  correlation  functions  and  the  cor¬ 
responding  potentials  of  mean  force  for  both  Coulomb  and  shielded 
fields. ^The  calculations  are  made  up  to  second  order  in  complexity. 
Corrected  expressions  for  the  potentials  of  mean  force  are  compared 
with  the  corresponding  Debye-Hiickel  results.  Significant  differences 
are  found;  this  is  particularly  true  for  the  shielded- field  cases. 

It  is  demonstrated  that  use  of  the  corrected  pair-correlation  func¬ 
tions  in  the  Moxer-Baranger  microfield  calculation  could  signifi¬ 
cantly  altar  the  results.  The  possibility  of  including  such  corrected 
results  in  the  line -broadening  theory  is  considered^  it  is  indicated 
that  the  influence  of  these  corrections  on  the  profi^Sif  the  Lyman-  yS 
line  should  be  experimentally  observable. 
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I.  INTRODUCTION 

The  theory  of  spectral  line  broadening  has  recently^  ^  been 
improved  to  the  point  where  the  profiles  of  both  overlapping  and  isolated 
spectral  lines,  arising  from  hydrogen  and  helium  plasmas,  are  well  under¬ 
stood.  Two  features  of  these  spectral  lines  are  frequently  of  special 
interest:  the  line  half  width  and  the  line  shift.  The  first  of  these 
quantities,  together  with  the  entire  line  profile  is  theoretically  pre¬ 
dicted  to  within  10%  of  the  experimentally  determined  values.  However, 
the  spectral  line  shift  is  generally  not  in  agreement  with  experiment; 
deviations  of  better  than  30%  are  frequently  found. Furthermore, 
these  variations  are  not  always  in  the  same  direction.  It  therefore 
remains  to  explain  these  shifts  satisfactorily. 

In  the  construction  of  the  general  line-broadening  theory,  numerous 
approximations  have  been  made.  A  review  of  the  underlying  assumptions 
leads  to  the  conclusion  that  the  most  probable  causes  of  the  shift  dis¬ 
crepancy  are:  (1)  the  crude  treatment  of  strong  electron  collisions, 

(2)  the  use  of  the  Debye-Huckel  pair-correlation  function  in  the  calcu¬ 
lation  of  the  electric  microfield  distribution,  (3)  the  approximate  treat¬ 
ment  of  ion  perturbations  in  the  cases  where  the  present  theory  is  not 
adequate. Although  these  quantities  will  be  discussed  in  detail  in 
succeeding  chapters,  it  may  be  mentioned  here  that  the  pair  distribution 
function  specifies  the  spatial  distribution  of  two  interacting  particles, 
while  the  electric  alcrofield  distribution  function  determines  the  spatial 
characteristics  of  the  electric  field  in  a  plasma.  Thus,  the  electric 
field  perturbs  the  radiating  atom  and  therefore  influences  the  nature  of 
the  radiation  itself.  It  is  known  that  in  some  instances  the  line  shift 
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is  particularly  sensitive  to  the  electric  field  distribution;  thus,  a 
method  of  checking  these  distributions  is  available. 
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Current  theories  of  line  broadening  rely  on  the  Debye-Huckel  distri¬ 

ct) 

bution  function  for  both  the  Coulomb  and  shielded  fields.  However,  as 
will  be  indicated  later,  the  assumptions  necessary  in  the  derivation  of 
this  function  are  not  applicable  to  many  plasmas  of  interest.  It  is  the 
purpose  of  this  thesis  to  calculate  the  pair  distribution  function  for 
the  Coulomb  and  shielded  electric  fields,  to  consider  in  detail  the  various 
orders  of  correction,  and  finally  to  indicate  what  significance,  if  any, 
the  deviations  from  the  Debye-Huckel  results  have  as  related  to  the  shift 
of  spectral  lines. 

For  the  purpose  of  this  investigation,  it  is  necessary  to  specify  a 
particular  plasma  model.  The  one  chosen  eliminates  infinitely  attractive 
potentials  from  consideration.  The  electrons  are  treated  as  negatively 
charged  particles  in  a  uniformly  neutralising  background,  and  the  ions  are 
assumed  to  be  shielded  by  electron  clouds  which  travel  permanently  with 
them.  That  such  a  model  is  applicable  to  practical  plasma  problems  may 
not  be  immediately  apparent;  further  discussion  and  justification  will  be 
presented  in  Chapter  III. 

The  motivation  for  choosing  this  particular  model  stems  from  the  work 
of  Moser  and  Baranger^  ,  which  describes  a  clever  method  for  calculating 
the  electric  microfield.  The  Baranger  approach  utilises  a  cluster-expansion 
program  which  is  reminiscent  of  the  well  known  Mayer  cluster  development. 
Particle  correlations  are  included  in  this  theory  as  a  perturbation 
expansion:  the  first  group  of  tens  contain  the  1-body  or  independent 
correlation  functions,  the  second  group  the  2-body  or  pair-correlation 
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functions,  etc.  For  such  a  process  to  be  valid,  increasing  orders  of 
correlation  must  become  very  small,  in  a  perturbation  sense.  Although 
both  the  electrons  and  ions  contribute  to  the  electric  microfield,  they 
do  so  in  a  very  different  manner;  in  view  of  this  difference,  discussed 
in  Chapter  III,  the  two  contributions  are  handled  separately.  The  relation 
of  the  corrected  correlation  functions  to  the  Mozer-Baranger  theory  will 
be  considered  later. 

In  order  to  effect  the  desired  correction  to  the  pair  distribution 

(2) 

function,  the  method  proposed  by  Salpeter  is  used.  This  procedure 

(52) 

has  previously  been  applied  by  Bowers '  to  calculate  the  Coulomb  poten¬ 
tial  of  the  mean  force;  however,  since  we  consider  several  aspects  of 
this  work  to  be  incorrect,  the  present  work  reverts  to  (2)  as  primary 
reference.  The  Salpeter  technique  embodies  one  of  the  cluster-expansion 
methods  frequently  used  in  statistical  mechanics.  It  differs  from  most 
other  approaches  in  that  it  Is  designed  for  use  with  systems  of  particles 
having  long-range  Interaction  potentials.  Although  the  most  notable  of 
these  is  the  Coulomb  potential,  the  method  should  also  be  useful  when 
dealing  with  potentials  which  are  long-range,  but  finite.  The  applicability 
criterion  is  that  the  potentials  involved  should  have  ranges  equal  to 

or  greater  in  magnitude  than  the  average  particle  separation,  given  by 
-1/3 

rQ  ~  (  p  -  particle  density).  In  the  cases  to  be  con¬ 

sidered  here,  this  requirement  is  satisfied;  it  is  automatically  satisfied 
for  the  Coulomb  case,  and  for  the  shielded  potential  it  is  also  valid  in 
the  range  of  densities  and  temperatures  of  interest  to  most  plasma  line 
broadening  applications.  Another  aspect  of  this  proposed  correction 
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procedure  is  Chat  it  involves  corrections  which  refer  to  particle  con¬ 
figurations:  The  first  term  of  the  pair-correlation  function  considers 
the  most  simple  clusters  of  interacting  particles;  the  second  correction 
takes  the  next  most  simple  configurations  and  so  on.  The  details  of 
this  expansion  technique  are  dealt  with  in  Chapter  IV. 

The  corrected  distribution  function  is  studied  as  a  function  of 
particle  density  and  temperature.  It  is  shown  that  in  the  given  range  of 
parameters,  the  second  order  (complexity)  contribution  is  much  less  than 
the  first  order  term.  However,  the  deviation  of  the  new  functions  from 
those  predicted  by  Debye-Huckel  theory  will  be  shown  to  be  quite  large 
in  some  cases.  In  light  of  this  last  disclosure,  the  expected  influence 
of  the  corrected  pair-correlations  functions  on  the  microfield  distri¬ 
bution  is  considered  and  the  need  for  a  revised  calculation  is  indicated. 

To  summarize:  In  Chapter  II  some  aspects  of  line-broadening  theory 
will  be  discussed  in  an  attempt  to  point  out  further  the  motivation  for 
the  problem.  Chapter  III  deals  with  the  electric  microfield  calculations 
of  Baranger  and  Moser  and  indicates  the  importance  of  the  pair-distri¬ 
bution  function  in  this  respect.  Salpeter's  technique  of  calculating  the 
pair -distribution  function  will  be  the  subject  of  Chapter  IV,  and  its 
application  to  both  the  shielded  and  Coulomb  cases  is  covered  in  Chapter  V. 
Chapter  VI  will  present  the  conclusions  of  this  study,  together  with 
recommendations  for  future  work. 


II.  SOME  ASPECTS  OF  LINE-BROADENING  THEORY 

From  1958  to  1962,  H.  &.  Griem,  A.  C.  Kolb,  and  M.  Baranger  developed 
a  theory  of  spectral  line  broadening  in  hydrogen  and  helium  plasmas.  This 
theory  pointed  out  the  importance  of  electron  broadening  and  provided  a 
coherent  treatment  of  the  joint  electron-ion  contribution  to  the  resultant 
spectral  line  profile. 

The  approach  is  simply  this:  electrons  and  ions,  treated  as  clas¬ 
sical  particles,  perturb  a  radiating  atom,  thereby  influencing  the  intensity 
distribution  of  the  resulting  radiation.  Looking  more  closely  at  the 
procedure  one  sees  that  in  general,  three  steps  are  followed: 

(1)  The  radiating  atom  is  assumed  to  lie  in  a  static  electric  field, 

E,  created  by  the  ions,  and  hence  is  subject  to  Stark  splitting. 

(2)  Electron  broadening  of  the  resultant  sharp  lines  is  treated 
by  the  impact  approximation. 

(3)  The  electron-broadened  spectrum  is  averaged  over  the  probability 
distribution  of  the  electric  field. 

Step  (1)  amounts  to  the  usual  Stark  effect,  where  the  electric  field  E 
is  calculated  for  a  particular  ion  configuration. 

In  carrying  out  step  (2),  the  electrons  are  treated  in  the  classical- 
path  approximation.  In  order  to  justify  this  assumption  that  the  electrons 
may  be  treated  as  classical  particles  with  straight-line  trajectories, 
one  must  state  that  the  momentum  transfer  in  such  an  impact  perturbation 

•■all*  !•«•*  1  where  Ap  momentum  transfer  and 

P 

P  =  initial  momentum.  Another  necessary  restriction  is  that  the  perturber 
wave-packet  dimension,  &  X,  be  small  compared  with  the  impact  parameter, 

X  Thus,  it  is  seen  that  we  are  dealing  with  interactions  which 
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on  the  average  are  weak  and  long  range. 

The  third  step  requires  that  an  electric  microfield  distribution  be 
calculated  for  the  plasma  model  previously  introduced. 

Fig.  1  illustrates:  (a)  the  nature  of  the  classical-path  approxi¬ 
mation,  and  (b)  the  type  of  overlapping  spectral  line  handled  by  this 
theory 


(a) 


(b) 


is  referred 


v 

to  as  the  collision  time. 


ii(t>  =  A+  v 


where  r. (t)  is  the  position 
vector  of  the  i  th  elec¬ 
tron  relative  to  the 
radiating  system. 


where  a  (  o<  ,  «',oc")  represents  an  initial  level  with  principal 
quantum  no.  "a"  and  substates  °<  ,  <*  ' ,  o<  "  and  b  (  j3  ,  ' )  represents 

the  corresponding  final  states. 


Figure  1. 
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The  transition  between  level  a(  a<  ,  oc  ' ,  o<  ")  ■■  ■—  V  b(  p  ,  *) 

may  give  rise  to  a  complicated  intensity  distribution  such  as  that  shown 
in  Fig.  2. 


J  <LJ) 


Figure  2.  Spectral  Line  Intensity  Distribution 

It  is  possible  to  derive  an  expression  for  the  spectral  intensity 
distribution  function  for  a  quan turn -mechanical  radiating  system  with 
time-dependent  dipole  moment,  ^  (t)  assuming  an  equilibrium 
situation. The  final  expression  for  this  quantity  is 


oO 


where 
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J  (6j)  —  Spectral  Intensity 
1a)  =  Angular  Frequency 
c  =  Speed  of  Light 
t  —  Time 

=  Boltzman  Factor  =  2  ^  e  ?  +  v(t)) 

2  =  Partition  Function 
k  1 

k  =  Boltzman  Constant 


T  =  Temperature  in  degrees  Kelvin 
Hq  =  Unperturbed  Hamiltonian 
V(t)  =  Time-Dependent  Perturbation 


ytl( t)  -  Dipole  moment  of  the  radiating  atom  at  time  (t) 
£  J  =  Average  over  all  perturber  configurations. 


The  average  specified  in  Eq.  (2.1)  refers  to  the  average  over  all  possible 
perturber  collisions,  e.g. ,  over  all  values  of:  (1)  impact  parameter; 

(2)  perturber  velocity;  (3)  collision  angles;  and  (4)  collision  time. 

This  averaging  process  is  necessary  because  ^ yu  (t) yU  (0)  is 
always  calculated  for  a  particular  instantaneous  perturber  configuration. 
Since  the  interaction  is  considered  weak,  it  is  assumed  in  the  following 
that  V  (t)  «  Hq  or  that 


p  i-  r1 P  % 


(2.2) 
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p0  is  independent  of  time  and  may  be  taken  outside  the  integral 


giving  (aside  from  a  constant  factor) 


The  problem  now  is  the  calculation  of  £  yU.  {t) ylA-  (0)  j  .  This  function 

is  frequently  referred  to  as  the  autocorrelation  function;  if  the  dipole 

moment  of  the  radiating  atom  is  /U  (t),  and  further,  if  the  perturbing 

/ 

interactions  serve  to  damp  the  dipole  moment,  there  exists  a  correlation 

time,  T  ,  past  which  the  time  integral  over  the  correlation  function 
c 

fast  approaches  zero.  The  correlation  time  then  sets  practical  limits 
on  the  time  integration  performed  in  the  evaluation  of  J  )  and  thus 
represents  a  measure  of  the  radiating  system  "memory".  In  general,  it 
is  not  possible  to  determine  this  function;  fortunately,  there  are  two 
special  cases  when  a  solution,  using  established  techniques,  can  be 
obtained.  These  are: 

(1)  The  collision  time  'Y  jy  T^  ,  the  correlation  time; 

(2)  The  collision  time  T  ,  the  correlation  time. 


The  first  case,  (1),  is  referred  to  as  the  quasi-static  approximation. 
Here  the  collision  is  so  strong  that  it  completely  disrupts  the  phase  of 
the  radiating  system  in  a  time  short  compared  to  the  collision  time,^  ; 
the  disruption  takes  place  so  rapidly  that  it  may  be  assumed  that  the 
perturbing  particles  (ions)  remain  fixed  during  \  .  This  first  approxi¬ 
mation  is  useful  for  dealing  with  ion  perturbations  over  major  portions 
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of  most  spectral  line  profiles.  It  Is  particularly  good  In  regions  well 
removed  from  the  center  of  the  line;  this  of  course  is  true  because  of 
its  strong-collision  nature. 

In  the  second  case,  the  impact  approximation  is  applicable;  here 
many  collisions  are  necessary  before  the  phase  of  the  radiating  system 
is  completely  disrupted.  This  implies  that  interactions  treated  in  this 
approximation  are  on  the  average  weak.  Electrons  are  very  successfully 
handled  with  this  technique  everywhere  but  on  the  far  wings  of  the  profiles; 
it  may  also  be  used  in  some  instances  to  determine  ion  contributions  at 
the  center  of  the  line. 

The  quasi-static  approach  is  used  to  accomplish  step  (1)  and  the 
impact  approximation  is  used  in  step  (2).  Hence,  at  the  end  of  these 
first  two  steps,  we  have  an  expression  for  the  spectral  intensity  which 
is  a  function  of  both  the  frequency  (a)  and  a  fixed  value  of  the  electric 
microfield,  B: 


J  ((a),  E)  =  [w2  +  (W  -d  -  C(B))2]  _1  (2.4) 

where:  w  =  line  half-width  as  determined  from  step  (2). 

d  -  line  shift  as  determined  from  step  (2). 

C(E)  =  a  function  of  the  field  strength  B;  it  is  proportional 

2 

to  B  for  the  linear  Stark  effect  and  to  B  for  the  quadratic 
Stark  effect. 

In  order  to  find  the  intensity  distribution  J  ((«>),  it  is  necessary  to 
average  J  ( u)  ,  B)  over  a  suitable  electric  microfield  distribution,  H  (E). 
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Therefore, 


oO 


or, 


(2.5  ) 


(2.6) 


Recall  Chat  Che  function  W  (E)  gives  the  probability  of  finding  an  electric 
field  of  magnitude  E,  at  a  particular  spatial  point  which  we  choose  to  be 
the  radiating  system.  This  function  is  strongly  dependent  upon  the  specific 
characteristics  of  the  particle  correlations.  The  exact  nature  of  the 
dependence  is  shown  in  the  next  chapter. 


III.  THE  ELECTRIC  MICROFIELD  DISTRIBUTION  FUNCTION 


A.  Introduction 

The  last  chapter  indicated  that  when  the  quasi-static  approximation 
was  applicable,  the  perturbers  involved  may  be  considered  to  produce  a  time- 
independent  electric  field.  To  arrive  at  an  electric  microfield  distribution 
function,  one  must  average  over  all  possible  perturber  configurations  or, 
what  amounts  to  the  same  thing,  over  the  possible  electric  fields. 

/A  \ 

Koltsmatk  first  calculated  this  distribution  function  for  a 
gas  of  charged  particles;  however,  in  doing  so  he  neglected  to  account  for 
the  fact  that  charged  particles  interact  with  one  another  and  hence  are 
strongly  correlated.  He  determined  the  probability  of  finding  a  particular 
electric  field  E,  due  to  a  static  random  distribution  of  singly  charged 

ions,  in  a  specific  range,  E  <  E  +  d  E.  Such  a  development  repre- 

o  o 

senta  a  strong-ccllis.ion  theory,  since  it  neglects  the  weaker  but  highly 
probable  perturbative  interactions  where  correlations  are  important.  There¬ 
fore,  it  is  not  surprising  that  the  Holtsmsrk  theory  incorrectly  describes 
the  center  of  spectral  lines,  but  gives  much  better  results  when  applied 
to  line  wings.  This  theory  will  also  be  briefly  mentioned  in  Chapter  IV. 

Since  this  initial  work,  others  have  attempted  to  improve  on 
the  theory  of  electric  microfield  distributions;  the  primary  effort  has 
been  to  correctly  Include  particle  correlation.  Ecfcer^\  and  Hoffman 
and  Theimer^®\  followed  the  Holtsmark  procedure  of  assuming  statistical 
independence,  but  replaced  the  Coulomb  field  with  a  shielded  field.  Broyles 
uses  the  Bohm-Plnes  collective-coordinate  technique  by  considering  short 
and  long  range  Interactions  separately  in  his  evaluation  of  the  microfleld 
distribution  at  a  radiating  ion.  Margenau  and  Lewis also  consider  the 
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electric  field  at  an  ion  due  to  perturbing  ions,  but  neglect  the  intension 
interaction.  This  last  method  also  has  the  disadvantage  of  not  being 
applicable  to  neutral  atoms.  Although  this  criticism  is  not  true  of  the 
Broyles®  approach;,  it  should  be  pointed  out  that  his  expression  for  the 
microfield  distributions  in  the  limit  X  approaches  infinity,  does  not 
reduce  to  the  expected  Holtsmark  result. 

In  two  papers,  ^Mozer  and  Baranger  put  forward  a  new  method 
of  including  correlations  in  the  calculation  of  the  electric  microfield 
distribution  functions.  The  construction  of  this  theory  involves  a  per¬ 
turbation  expansion  and  has  the  desirable  feature  of  reducing  to  the  Holtsmark 
result  in  the  limit  where  T  approaches  infinity.  To  be  practical,  we  must 
restrict  its  use  to  plasmas  where  the  interaction  energies  are  considerably 
less  than  the  thermal  energies,  i.e. ,  less  than  kT.  It  is  possible  to 
calculate  all  orders  of  correlation,  but  only  the  two-body  example  is 
developed  in  detail. 

The  actual  calculation  of  the  miorofield  distribution  function 
will  Involve  the  Fourier  transform  of  the  distribution  function  rather  than 
the  function  itself;  here,  as  in  many  instances,  the  transform  is  easier 
to  handle.  In  the  following  work,  E,  refers  to  an  arbitrary  electric  field; 
it  will  later  be  used  to  refer  to  both  the  Coulomb  and  shielded  fields. 
Following  tbe  notation  of  Mozer  and  Baranger,  we  refer  to  the  collective 
electron  field  as  the  high-frequency  field,  and  the  resultant  ion  field  is 
termed  the  low-frequency  field. 

B.  The  Mozer.  Baranger  Theory^1" 

Consider  a  volume  V  with  a  given  number  of  particles  K,  the 

.,  r„ . .  r  ;  the  particle  density  is  given 

~2  ~n 


coordinates  of  which  are  r 
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by  P  =  S  .  Furthermore,  It  is  assumed  that  thermal  equilibrium  conditions 
prevail.  The  quantity  to  be  calculated,  W  (E),  gives  the  probability  o£ 
finding  a  total  electric  field  E  in  a  range  E  to  E,+  dE_at  a  given  point 
in  space.  It  is  usually  convenient  to  take  this  point  as  the  origin  of 
the  coordinate  system.  The  Fourier  transform  of  the  field  distribution 
function  is 


F  (k)  = 


ik.g. 

e  ~  W  (E)  dJ  E 


(3.1  ) 


where  the  field  E  at  the  origin  is  equal  to 

N 

*•  -  i u 

i=i 


E-  =  field  produced  by  one  of  the  particles  located  at  point  r^.  Mow 
It  is  well  known  that  Eq.  (3.1)  may  also  be  expressed  in  terms  of  the 

particle  distribution  function  P  (j^f  •••  r^,).  Thus,  we  replace 

3  o 

W  (E)  d  E  in  the  Fourier  transform  by  P  (r,  • • •  dJr-  where 
^  ''l  H!'  M 

P  (jr^  •••  r^  )  gives  the  probability  of  finding  the  M  particles  in  a 
particular  configuration.  For  a  classical  gas  in  thermal  equilibrium, 

P  (£]_  •••  £d)  Is  the  Boltzmann  distribution. 


P  (r 


■"  V 


_  e 


/'**/e=/SU  *’*  V  -jj-  d3  rt 

1=1 


volume 


(3.2) 
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Here,  we  have  defined 


.3  .3 

dr,  • • •  drJ 
1  N 


With  the  above  substitution,  Eq.  (3.1)  becomes 


%«!+—  +  Sn)  JL  , 

F(k)  =  •••  e  p(£l  •**  JEjj)  '*  d  ri  (3.3) 

i=l 


Using  a  trick  similar  to  that  employed  in  the  usual  cluster-expansion 
procedure,  we  may  write: 


“  ‘  Ei  % 

e  =  e  -1  +  1 


(3.4) 


=  Pt  +  1 


Substitution  of  this  expression  into  equation  (3.3)  gives: 

KW  -[■■■[  [i  +  ?  P  +  fi  $  *  -]  (3.; 


- »  <y +  ?!  (y +  ••• + '  <y +  ••• 


(3.5a  ) 


where 


(k) i <^rr< F(ii- •  -a>Y > 


Since  P  (r^  -  • •  r^)  £s  symmetric  with  respect  to  an  interchange  of  particle 

coordinates,  it  is  permissible  to  write  F  (k)  as 

m  ~ 


•5»  >  a3*. 


(3.6) 


Note  that  in  this  expression  it  is  possible  to  integrate  over  all  but  the 
specified  M  particle  coordinates  thus  obtaining  a  new  probability  function 
which  depends  only  on  the  M  coordinates  (r^  •••  r^) ,  and  is  equal  to 
VM  times  the  M-body  correlation  function.  Hence,  we  write: 


FM<i9 


N  ! 

M!  (N  -  M)! 


(3-7) 


P  can  now  be  expanded  in  a  perturbation  series  involving  increasing  orders 

n 

of  correlation.^®)  Thus,  the  first  term  contains  only  theproduct  of  one- 
body  or  Independent  probability  functions;  the  second  term  involves  the 
two-body  Interaction  corrections,  while  each  successive  term  includes  cor¬ 
rections  involving  one  more  particle  added  to  the  particle  cluster.  To 


illustrate: 
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where  s^  (£^)  represents  the  one-body  or  independent  probability  function, 
is  the  two-body  correlation  function,  and  so  on.  Since  the  s  func¬ 
tions  are  also  symmetric  with  respect  to  an  interchange  of  particle 


coordinates,  F  (k)  becomes 
M 


F  (k)  =  - NJ - 

M! (N  -  M)! 


i=l 


u  i  trH  I  C  C  *1 

^ - —  <C(M)  JTCp  8  (r  )  dr 

M!  (H-H)!  j  1  ...  1  11  '  i  i  i 

J  J 


+  C(M) 


FL  77  *1  <ij>  'a'n-  tU  d\ 


n 

:1<M,.1,  2,  J-Jv  W-y 


4*  •  •  • 


+  C(M) 


J-Jtt  *!%>  Er  Ej'TPPi 


H 
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The  coefficients  C  (MV  are  numbers  equaling  the  number  of 

ways  to  choose  distinct  groups  of  y  ' '  Particles  out  of  a 

total  of  M  identical  particles,  with  <n  ”  *  ^ •  Distinct 

groups  implies  that  all  the  particles  of  a  specific  group  are  different 

from  all  particles  of  any  other  group.  In  general,  if  there  are 

^i  J  etc*>  with  the  condition  that  >iA+M  +  " 

+  X  M  =  M,  the  coefficient  is 
r  /  r 


C(M) 

A  ”■  A’  **  j/r***  A»  •' 


(/^Z  1  ...  (3.10) 


One  can  define  a  quantity  h  (k)  as 


<y  mf-fpi<r2  -  Vi  <c  * 


(3.11) 


With  this  definition,  the  Fourier  transform  of  the  field  distribution, 
in  the  limit  as  N  and  V  approach  infinity,  while  the  density  P -  £ 
remains  constant,  is 

°o  \  \ 

F(k>  =  [  2  C(M)  h  00  **•  h  jk)l 

all  partitions  r 
of  M 


-f  Vl(y]V 

1  \=0  A  :  Al  'd  •>'1! 
/* 

OO  r  /jl 

v-id 


(3.12) 
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=  exp 


This  expression  for  F  (k)  is  one  frequently  found  in  the  conventional 
cluster  expans ions (^0) ^  an  exponential  of  a  power  series  where  successive 
terms  contain  higher-order  clusters.  It  can  be  shown  that  in  the  high- 
temperature  limit  only  h^  (k)  remains  and  that 


which  with  E_(r)  -  —  can  be  shown^D  to  be  the  Fourier  transform  of 

r3 

the  Holtsmark  distribution  function. 

Now  consider  the  ptb  term  in  the  series  argument  of  the  exponential, 

p 

hp  (k)  ;  this  may  be  expressed  as  an  Integral,  as  indicated  by 

Eq.  (3.11),  the  value  of  which  is  determined  by  the  properties  of  (p ' 

and  Sp  .  The  salient  features  of  the  sp  are:  (1)  They  are  different 

from  zero  only  when  the  p  particles  are  clustered;  (2)  The  region  of 

particle  correlation  is  finite  due  to  screening  by  the  charged  particles 

in  the  gas;  (3)  sp  goes  to  zero  as  T  approaches  infinity.  Since  the  (p ^ 

decrease  rapidly  as  the  distance  from  the  field  evaluation  point  is  increased, 

the  integral  over  all  space  is  convergent.  For  finite  values  of  k,  hp(k)  is 

also  finite.  That  the  convergence  of  the  integral  is  dependent  on  the 

factor  is  apparent  since  for  large  values  of  all  the  r,  the  s  does  not 

~i  p 

necessarily  vanish  since  the  particles  may  still  be  close  together.  One 


would  also  expect  that  for  finite  k  and  large  temperatures  the  successive 
terms  in  the  exponential  of  (3.12)  would  be  proportional  to  increasing 
powers  of  the  inverse  temperature,  then  only  a  few  terms  in  the  expansion 
are  necessary  to  get  a  good  approximation  for  P  (k)  in  the  high-temperature 
region. 

Eq.  (3.12)  is  quite  general  and  may  be  of  interest,  in  the 
high-temperature  limit,  for  a  number  of  fields.  The  electric-field 
distribution  follows  immediately  upon  applying  the  inverse  Fourier  Trans¬ 
form  to  (3.12), 

.  C  -ik  •  E 

W(E)  =  -i —  |  e  ~  ~F(k)  d3k  3.13 

~  (2TT)3  J 

C.  The  High-Frequency  and  Low-Frequency  Pields 

Before  continuing  in  the  calculation  of  the  electric  field  dis¬ 
tribution  function,  we  will  discuss  the  physical  nature  of  the  electric 
field  contributions.  Due  to  the  differences  in  the  time  variation  of 
their  resultant  fields,  it  is  desirable  to  treat  separately  the  electronic 
and  ionic  contributions  Do  the  total  field.  He  follow  the  notation  of 
Mozer  and  Baranger,  referring  to  the  field  of  the  electron  as  the  high- 
frequency  component  and  that  of  the  ions  as  the  low-frequency  component. 
That  it  is  reasonable  to  consider  the  contributions  separately  is  indicated 
not  only  by  the  differences  in  the  time  characteristics  of  the  fields, 
but  also  by  the  fact  that  the  collective  fields  of  the  electrons  and  ions 
are  much  different.  In  order  to  discuss  these  matters  further,  it  is 
necessary  to  establish  suitable  time  scales  which  are  characteristic  of 
the  ionic  and  electronic  motion.  For  this  purpose,  the  concept  of 
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relaxation  times  are  introduced.  If  we  consider  a  gaseous  system  of  like 
particles  which  originally  have  an  arbitrary  velocity  distribution,  the 
time  required  for  a  Maxwellian  velocity  distribution  to  be  established 
through  collisions  is  referred  to  as  the  relaxation  time. 

It  has  been  shown (53)  that  in  a  positive  ion-electron  gas, 
where  the  velocity  distributions  of  both  electrons  and  ions  are  initially 
arbitrary,  the  electrons  will  come  to  Maxwellian  equilibrium  much  more 
rapidly  than  will  the  ions.  The  appropriate  relaxation  times  which 
govern  this  approach  to  the  Maxwellian  distribution  are  taken  as  charac¬ 
teristic  of  the  ionic  and  electronic  motion. 

Returning  now  to  a  consideration  of  the  physical  nature  of  the 
two  field  components,  we  realize  that  it  is  not  correct  to  use  a 
Coulomb  field  for  the  individual  ions;  such  a  procedure  would  completely 
neglect  the  ion-electron  interaction.  The  low-frequency  component  contains 
a  contribution  from  the  electrons  since  the  ion-electron  Interaction 
establishes  a  shielding  electron  cloud  about  each  of  the  ions.  Thus, 
the  low-frequency  component  consists  of  a  sum  of  all  the  separate  shielded 
fields.  The  correct  method  for  obtaining  the  low-frequency  component 
is  to  calculate  the  time  average  of  the  electric  field  for  a  given  ion 
configuration;  the  time  over  which  the  average  is  computed  should  be  short 
compared  to  the  ionic  relaxation  time,  but  long  compared  to  the  cor¬ 
responding  electron  time.  Since  we  are  involved  with  a  time-averaged 
component  here,  it  is  proper  to  use  a  Debye-Huckel  shielded  field  for 
each  ion. ^  Note  that  the  shielding  is  due  solely  to  electrons;  ion- 
ion  interactions  are  taken  into  account  in  the  pair-correlation  functions. 
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Because  the  ionic  time  of  relaxation  is  much  longer  than  that 
for  electrons,  an  average  of  the  high-frequency  component  of  the  micro¬ 
field  over  times  the  order  of  the  ionic  relaxation  time,  will  vanish. 

Hence,  it  is  reasonable  to  consider  the  high-frequency  component  as  a 
gas  of  electrons  imbedded  in  a  uniform  neutralising  background. 

The  above  procedure  is  essentially  that  followed  by  those  who 
describe  plasma  behavior  by  collective  coordinates.  However,  the 

arguments  presented  here  are  not  intended  to  imply  anything  more  than 
plausibility.  The  concept  of  low-frequency  and  high-frequency  com¬ 
ponents  can  be  derived  more  rigorously  from  the  collective-coordinate 
approach. 

At  this  juncture  the  Mozer-Baranger  approach  resorts  to  the  usual 
Debye-Htickel  method  to  determine  the  shielded  fields  of  the  individual 
ions.  The  validity  condition  for  this  technique  is  that  the  volume  of 

the  electron  cloud  is  large  compared  to  the  average  volume  per  particle. 
Thus,  the  ion  field  is  found  to  be 


(3.14) 


where  A  =  Debye  length. 

Consider  now  the  two-body  correlation  function  which  will  sub¬ 
sequently  be  used.  For  the  high-frequency  component,  treated  as  a  gas  of 
electrons  against  a  uniform  neutralising  background, the  potential  of  the 
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and  the  corresponding  correlation  function  is: 


i2)  = 


Z3  W 


5a) 


(3.15b) 


These  expressions  may  be  derived  by  the  Debye -Huckel  method,  but  for 
the  purpose  of  this  dissertation  they  are  defined  and  derived  from  a 
different  point  of  view  in  the  next  chapter,  as  are  the  following  results 
for  the  low-frequency  component: 


2  -  i!i- 

w  (r  ,  r  )  =  £ -  e  AVT 

2  M 

ri2 


and 


I2>  *  * 


(3.16a) 


(3.16b) 


It  la  seen  from  these  equations  that  the  only  difference  between  the 
two-body  correlation  functions  for  the  electron  and  ion  cases  is  contained 
in  the  definition  of  the  shielding  length,  X  or  X  /{ 2  '. 

D.  The  Fourier  Transform  of  the  Field  at  a  Neutral  Point 

The  following  procedure  will  be  considered  in  some  detail  and 
is  applicable  to  either  the  hlgh-or  low-frequency  case. 

The  neutral  point  at  which  the  field  is  to  be  evaluated  (the 
field  evaluation  point)  is  chosen  to  be  the  origin  of  the  coordinate 
system.  In  this  case  the  s^  (r^>  =  1  and  the  one-body  term,  ^h^(k)  ,  is 

f  ik.E  , 

P\  qo  -fi  J  c  ~  n  - 1)  A, 


(3.17) 
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Since  E^  =  E^(r^)  we  “ay  set  k  =  k  with  no  loss  of  generality, 
and  Eq.  (3.17)  reduces  to 


ao 


y^h^(k) 


47rr  (1  "  J0(kEl))  ri  % 


(3.18) 


where  j  is  the  zero-order  spherical  Bessel  function and  E.  is 
the  magnitude  of  the  single-particle  field;  a  shielded  field  is  used  for 
the  low-frequency  case  and  a  Coulomb  field  for  the  high-frequency  case. 
In  Idie  latter  instances  (3.18)  may  be  evaluated  analytically,  yielding 


/?hl(k)  ='^P(27r  ek>' 


(3.19) 


which  depends  only  on  Che  magnitude  of  k.  Prior  to  writing  a  similar 
expression  for  the  low-frequency  case  (L.F.)  it  is  convenient  to  make 
the  following  definitions: 

(a)  A  unit  of  length  r  is  defined  by 

o 

.  3/2  , 

A  <2ir> 

i 

r  is  almost  equal  to  the  radius  r  ,  where 
o  o 

-  * 

(b)  The  "normal"  unit  field  strength  E  is  defined  by 

E 


o 


.e. 

.2 

o 
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(c)  Two  additional  units  are 


$  = 

E 


=  J  k  |  E 


These  units  will  be  employed  in  the  remainder  of  this  chapter.  Thus, 
(3.19)  becomes: 


yOhj^  (k)  =  -  x 


(3.19a) 


The  one-body  term  for  the  shielded  field  is  given  by 


O  h  (k)  =  -  - 

7  1  (8fr)l/2 


oo 

j  b + 


Z2  d  z 


(3.20) 


where 


c  -2  1/2  T  r  1/2 

|  =  Z  (1  +  ro/A  )  x  Z  exp 


)  x  Z 


Z  =  r  coordinate  divided  by  (ek) 


Further  let  us  define 


_  ro  1/2 

y~  x x 


(3.21) 


Hence,  Bq.  (3.20)  may  be  written  as 


(k)  =  -  x  \f>l  (y) 
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This  relation,  which  defines  (y),  must  be  evaluated  numerically. 

Now  the  two-body  term,  h  ,  is  found  by  letting  U  =  2  in 

Z  '  i 


Eq.  (3.12).  Thus,  we  have 


h2<k)  -  . 


—  ^2  s2<ll -ij)  "S  “3 


r2  (3.22) 


where  from  Eqs.  (3.8)  and  (3.15) 


*2  (V  V  = 


(3.23) 


In  the  high-temperature  limit  one  may  linearize  this  expression: 


(3.24) 


which  implies  that 


/>2  -  /°2 

S'  h,  (k)  =  _c_ 

2  2  2!  kT 


?1^2  w2  <£i»  i2)  ^  d3r2  5  O- 25a) 


or  in  dimensionless  units. 


Pl  ~  3/2 

l  h2  =  x  ^2 


(3.25b) 


Here  (^(y)  contains  the  two-body  correction  term.  Thus,  jp  (y)  and 
jp  2  (y)  refer  to  the  appropriate  one-body  term  and  the  two-body  cor¬ 
rection  respectively.  For  the  very  high- temperature  cases  where  the 
function  w2(rj,  r^)  is  given  by  the  Debye-Huckel  theory,  it  is  possible 
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to  evaluate  this  expression  by  expanding  the  various  terms  in  the  inte- 
grand  in  spherical  harmonics.  However  we  are  not  interested  in  this 
solution  here.  Attention  will  be  confined  to  comments  about  the  general 
character  of  Eq.  (3.24)  and  its  application  to  the  calculation  of  micro¬ 
field  distribution  function.  It  has  been  shown  that  the  value  of 
(k)  for  the  high-frequency  component  is  much  larger  than  the  cor¬ 
responding  quantity  for  the  low-frequency  case;  this  is  of  course  due 
to  shielding  effects. 

The  Fourier  transform  of  the  field  distribution  at  a  neutral 
point,  Eq.  (3.12),  to  the  order  of  approximation  used  here  can  be 
expressed  in  dimensionless  units  as 


F  (x)  =  exp 


(y>  -  f2  (y> 


(3.26) 


An  example  of  the  two-body  correction,  y  (y)  for  the  low-frequency 
case  is  shown  together  with  the  corresponding  one-body  term  for  com¬ 
parison  in  Fig.  3.  From  this  graphic  illustration  it  is  apparent  that 
for  values  of  x  which  do  not  make  the  exponential  exceedingly  small, 
and  for  values  of  rQ/  X  from  zero  to  one,  ^  ^  (y)  indeed  much 
smaller  than  ^  (y).  When  the  two-body  correction  is  small  compared 
to  the  one-body  term,  corrections  from  higher-order  correlations  should 
be  negligible  for  the  range  of  parameters  in  which  we  are  interested. 

E.  Microfield  Calculation 

The  Fourier  transform  of  the  field  distribution,  F  (k),  was 
shown  to  be  independent  of  the  k  direction.  This  implies  that  the  dis¬ 
tribution  function  itself  should  also  be  independent  of  the  k  direction; 
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Several  graphs  representing  H(£)  for  various  cases  and  values 

of  r  /  \  have  been  prepared.  However,  only  the  graph  In  Fig.  4  will  be 
o' 

discussed.  Here,  for  rQ/ X  ~  0.6,  the  low  frequency  H(£)  is  shown 
for  two  cases:  (1)  only  one-body  terms  included,  (2)  one-  and  two-body 
terms  considered.  Note  that  even  though  ^  (y  ) »  its 

inclusion  has  a  pronounced  effect  on  H(S),  the  microfield  distribution 
function;  specifically,  it  raises  the  distribution  maximum  and  shifts 
it  markedly  towards  weaker  fields.  The  conclusion  is  that  even  though 
\j)  2  (y)  is  *  very  small  quantity,  its  nature  strongly  influences  the 
form  of  H(&).  It  will  be  shown  in  chapter  V  that  the  corrected 
is  considerably  different  from  the  Debye-Huckel  result.  Hence,  this 
correction  might  be  expected  to  alter  W  (E)  or  H  (o)  considerably. 


* 

1 
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IV.  CALCULATION  OF  PAIR-CORRELATION  FUNCTIONS 
A.  Introduction 

A  large  number  of  problems  in  statistical  mechanics  is  con¬ 
cerned  with  calculation  of  correlation  functions  for  systems  with 
large  numbers  of  particles.  The  method  of  cluster  expansions (23) 

has  been  most  successful  in  this  regard  when  applied  to  systems  of 
particles  with  very  short  range  interactions.  Using  this  technique, 
one  is  able  to  derive  infinite  expansions  in  powers  of  the  particle 
density,^-5  ,  for  the  many-body  correlation  functions.  If  the  density 
of  the  system  is  sufficiently  low  and  if  the  particle  interactions  have 
short  range  these  expansions  converge  rapidly;  this  is  true  when  dealing 
with  neutral  molecules. 

However,  in  the  case  of  long-range  potentials,  of  which  the 
Coulomb  interaction  is  the  prime  example,  this  cluster  expansion  approach 
is  not  applicable  regardless  of  how  dilute  the  system  may  be.  Debye 
and  Hu'ckel  proposed  a  method  (1923)  for  treating  very  low  density  systems 
involving  Coulomb  interactions.  Since  this  work,  several  authors (22) (47 ) 
have  developed  techniques  for  use  with  systems  of  finite  density  which 
yield  for  the  Debye-Huckel  result  as  a  first  approximation  and  which, 
at  least  in  principle,  also  offer  higher-order  correction  terms. 

Generally,  these  improved  methods  express  the  correlation  functions  in 
terms  of  complicated  integral  equations. 

The  present  chapter  will  essentially  present  the  Salpeter 
approach  to  the  problem  of  calculating  correlation  functions  for  a 
collection  of  particles  which  interact  through  long-range  potentials. 
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The  problem  of  how  to  regroup  the  terms  of  a  cluster  expansion  so  as  to 
give  meaningful  results  when  applied  to  situations  involving  long  range 
potentials^  is  considered  in  detail.  In  this  development  of  the 
expansion  approach  it  will  be  assumed  that  we  are  dealing  with  an 
equilibrium  system  of  volume  V  which  contains  N  particles  of  mass  M; 
there  exists  a  central  interaction  potential  U(r^)  acting  between 
any  pair  (i,  j)  of  particles.  The  theory  is  classical  in  nature  and 
assumes  only  two-body  interactions.  In  order  that  the  classical  assump¬ 
tion  be  correct,  infinitely  attractive  potentials  are  forbidden. 

B.  Preliminary  Statistical  Mechanical  Considerations 

In  this  section,  some  general  concepts  from  statistical 
mechanics,  which  will  be  necessary  in  the  discussion  of  the  Salpeter 
theory,  are  discussed.  Of  central  importance  in  statistical  mechanics 
is  the  partition  function.  It  is  defined  as 


Z  = 


h3*^! 


_  _H_ 

e  kT  ^3  .3  3 

*  (in.  ...  d  n 


Pl  •••  d  pN  d  rx  •••  d3rM  (4.1) 


where: 


d  pt  =  d  p*t  dpyI  dp,i 


d  rt  =  d  d  dz1 


N 


2 

“  "  Y,  af  (p«l  +  pyl2  +  pzl2)  +  Hll  •••  V 

i=l 


H  is  the  Hamiltonian  of  the  system. 

U  is  the  potential  energy  of  the  system. 
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p.  =  momentum  of  the  ith  particle. 

=  coordinate  of  the  ith  particle, 
h  =  Planck's  Constant 


Note,  also,  that  N  and  V  are  considered  large;  it  is  frequently  con= 

venient  to  have  both  V  and  N  approach  infinity  in  such  a  way  that  the 

density,  P  =  H.  ,  remains  constant. 

'  V 

Explicitly  using  the  above  H  and  carrying  out  the  integration 
over  the  momentum,  we  have 


Z  = 


_  N 

Qn  v 

N  !  A3N 


(4.2) 


where 


%  =  V 


-N 


i 


<£n; 


d3r 


d3r„ 


h 

(2TTmkT)1/2  * 


u(5i  •••  Jn>  -  ~ 

kT 


U(r 


■Zl 


(4.3) 


Now  is  referred  to  as  the  configuration  integral.  Making  the 
assumption  that  the  potential  energy  u(r^  •••  r^)  ts  equal  to  the  sum  of 
pair  potentials, we  write 


■Si 


)  = 


1  ±  i<  j£ 


"(V 


.  zz 


"ij 


(4.4) 


i£.i<  j£.N 


Then 
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*N  = 


Uj  iJ 


d3r. 


dV 


(4.5) 


(1)  Distribution  and  Correlation  Functions^) 

It  is  now  possible  to  define  a  distribution  function. 

Suppose  we  observe  a  particular  equilibrium  configuration  of  N  particles 

whose  interaction  potential  is  given  by  u  (r,  •••  r  )„  The  probability 

~1  ~N 

3  i  3 

that  particle  1  will  be  in  d  r^,  2  in  d  r^,  and  N  in  d  rN>  is  given  by 


PN  <£i  •“  £M>dri***<irM 


-N  -1  _u(£1*  •  4E||) 

V  %  e  d3rL  ...  dJrN 


(4.6) 


The  probability  p  (r.  •*.  r  )  d3r.  •  •  *  d3r  that  particle  1  will  be  in 
n  i  n 

03  3 

dJr, ,  2  in  d  r„  n  in  d  r  regardless  of  the  arrangement  of  the 

L  *■  n 

remaining  N  -  n  particles,  is  the  sum  of  all  probabilities  (4.6)  con¬ 
sistent  with  the  specified  (1  .•*  n)  configuration:  i.e. , 


P  (r, 
n  '-1 


r  ) 


-N  -1 
v  % 


n+1 


(4.7) 


subject  to  the  normalization  condition 


d3r 


(4.8) 


If 

2  in  d3r  , 

2 


we  remove  the  restriction  that  particle  1  be  in 
etc. ,  and  specify  only  that  one  particle  be  in  each 
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element  dr  ,  the  correct  probability  function  to  use  is  P  (r,  • • •  r  ) 
i  J  n  m.  mi 


it  is  related  to  p  (r  •••  r  )  by 
rn  V-l  1 


P  <*1  •”  =  -J3-L 

7  n  ^  ^  (N-n) 


P  (*i  •••  jl) 


where  the  factor  - : -  gives  the  number  of  possible  ways  such  an 

(N-n)  ! 

arrangement  can  be  made.  The  normalization  condition  here  is 


j-  V-  «■ 


4>  iri 


3  3  m  » 

r  •  •  •  d  r  = 


(4.9a) 


(N-n)! 


Now,  since  this  dissertation  is  concerned  with  plasmas  , 
let  us  consider  the  application  of  these  distribution  functions  to 
gaseous,  or  more  generally,  fluid  systems.  As  examples,  we  choosey*3 ^ 
andy^2*  From  Bq.  (4.9) 


P\  (rl)  =  H  Pi  <rl> 


-  N  Qh_1  v"H  j  •••  e 

4 


<*1  ,  &  =  N  (N  -  1)  Pl  (r1,  £2) 


-1  -N 

=  N  (N  -  1)  0  V 


u(*l  .3 

dr  •••dr 


'u(£i  £n)  3 


(4.10) 


(4.11) 


The  normalization  condition,  Iq.  (4.9a),  with  n  =  1  yields 


J 7*1  ^  A>  =  * 


(4.12) 
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boundary  layer; 


In  a  fluid  all  points 
that  is,  (r^)  Is 


r  are  equivalent  if  we  neglect  the 
^•1 

independent  of  r  .  Hence,  from  Eq.(4.12) 


N 

V 


the  particle  density 


The  two-particle  distribution  function  for  the  fluid  case  is  only  a 
function  of  the  distance  between  two  points.  This  may  be  expressed  as 


Pi  ''ll  ■  lP  ■  Pi  <Sl2) 


Again,  using  Eq.  (4.9a)  we  find. 


r  r 


P 2  <ri  »  la}  d3r  d3r  =  - SJ - 

/  ~l  **  1  2  (H  -  2)  ! 


or 


ft  <‘u>  A*  - 


An  examination  of  Eq.  (4.8)  shows  that  if  the  distribution  of  particles 
is  completely  random, 


Pn  =  V 


-n 


(4.13) 


p  and  P  are  said  to  be  "independent”  if  they  satisfy  the  following 
n  /  n 


relations: 


pn  'll  -■  In>  "'1  (il>  •••  pl 


Pn  'l.l  ”•  ‘Pi  (Sl>  '"Pi  fc>> 


(4.14) 
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If,  however,  the  n  probabilities,  P.  (r  ),  are  not 

i 

"independent",  it  is  possible  to  define  a  correlation  function,  g^, 
which  provides  a  measure  of  the  deviation  of  p  or  P  from  the 
"independent"  value.  Considering  only  ^  ,  we  can  write 

f*  •••  Zj  7^  (£i>  —A  (Ji 


Returning  to  the  fluid  case, 


(4.15) 


n  n 

y°n  <Xl  '•*  In*  =  f>  8n  (£l  £nJ 


From  Eqs.  (4.7),  (4.9),  and  (4.15), 


*  u  ’  "I  -u(r  . ..r  )  ,  1 

«n  “•  A,)=  -  Si  •“  e  d  r  ••*d  rN  <4'16) 

»°  (H  -  n)!  *  ir+1 


Since 


U -  _  i  _  n(n-l)  +  .. 


Nn(N-n)! 


we  can  simplify  Eq.  (4.16)  when  terms  of  order  i  are  negligible. 

N 

Thus,  we  write 

n-N  -if  [  -u(r  •  •  mO  , 

*o  ‘il  V  *  v  %  "■  *  “  4  *.«  ■"  4  r»  (4'17) 


=  V  P  (£,  ...  r  ) 
n  n 
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Closely  connected  with  the  correlation  function  is  the  potential  of  the 

(3) 

mean  force,  defined  by  the  equation'  ' 


I  (r,  ...  r  )  =  exp  J  -  w  (r.  •••  r  )  [  , 


(4.17a) 


where  w  (fc,  •••  r  )  is  the  potential  of  the  mean  force  expressed  in 
n  ^1  /wn 

units  of  kT.  Now  recall  that  we  are  considering  a  system  of  N  particles 

which  interact  through  a  potential  u(r  ...  r  ),  and  that  of  these,  only 

~1  Ml 

n  have  their  coordinates  (r,  • • •  r  )  specified.  For  such  a  system  it 

(3) 

can  be  shown  that  the  appropriate  negative  gradient  of  w  (r  •  •  •  r  ) 

n  ^  a  n 

gives  the  mean  force  acting  on  any  one  of  the  n  particles,  e.g. ,  particle 
oC  .  The  mean  potential  differs  from  ordinary  static  potentials; 
in  addition  to  the  usual  contribution  from  the  n  specified  particles, 
it  also  has  a  contribution  from  the  remaining  (N  -  n).  Hence,  if  we 
hold  the  n  particles  fixed  and  average  over  the  other  (N  -  n),  we  find 
that  the  mean  force  acting  on  particle  oC  is 


(2)  Cluster-Integral  Expansions 

In  equation  (4.5)  there  is  a  term  in  the  Integrand, 
exp  [-!“»]•  Now  introducing  the  function  it  is 


possible  to  write: 


-2u 


u  =-rre 


=  TT(i  +  tLi) 

lii<  j£-N 


(4.18) 


i  +  Z.  +  2.Z  f,.  fkl  +  ••• 

ij  ij  ij  kl  kl 
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where  (i ,  j)  refers  Co  any  one  of  Che  possible  N  (N-l)/2  pairs.  Each 
term  In  this  expansion  can  be  represenCed  graphically  by  a  diagram 
which  corresponds  Co  a  facCor  f  (r^);  Che  CoCal  number  of  connecCing 
lines  in  Che  graph  equal  Che  number  of  facCors ,  f  ,  in  Che  parcicular 
term  under  consideraCion.  NoCe  ChaC  for  a  finite-range  potential  f 
goes  to  zero  rapidly  as  Che  separation  r  becomes  large.  To  illustrate 
the  cluster  principle  further,  consider  Che  various  possibilities  for  a 
set  of  three  particles: 


Figure  5  Clusters  of  Three  Particles 

The  terms  and  corresponding  diagrams  comprising  the 
expansions  (4.18)  can  be  classified  and  simplified  systematically.  For 
any  diagram  there  is  a  definite  number  of  points, ,  which  are  at  least 
singly  connected  to  particle  1.  The  portion  of  the  diagram  which  involves 
the  +  1  connected  points  is  designated  a  cluster;  the  remaining 
N  "  ($.+  1)  points  are  not  connected  with  the  given  (  X  +  1)  cluster, 
and  may  be  arranged  in  an  arbitrary  manner. 

How  consider  the  contribution  to  from  a  configuration 
involving  an  (  +  1)  particle  cluster  with  the  remaining  N  -  (  Jl  4-1) 

particles  arbitrarily  distributed. 


The  contribution  from  the  (1+1)  cluster  is  given  by 


IX+1 


-  <$+!> 


=  V 


dr  ...  d3r  Tf  f 

1  i+  1 


(4.19) 


where  the  product  of  f  includes  all  such  factors  occurring  in  the 

diagram.  ^  multiplied  by  an  integral  over  the  remaining  N  -  (  1+1) 

particles  yields  the  contribution  to  from  the  above  arrangement.  The 

sum  of  all  possible  distributions  of  the  N  -  (1+1)  particles,  while 

holding  the  (1+1)  fixed,  gives  a  total  contribution  to  Q  of 

N 


I  Q 

l+l  H  -  (1+  1) 


where 


N  -  (1+1) 


=  V 


-  H  +  U  +  1) 


d3r 


1  +  2 


3 

•  •  •  dr  e  1J 


(4.20) 


Multiplying  by  a  factor  -(ft  ~  O'  allows  for  the  fact  that  any 

ilD*  -  (l+i)]! 

of  the  (N-l)  particles  can  be  used  together  with  particle  1  in  making 


up 


the  (1  +  1)  cluster.  At  this  point,  we  define  a  cluster  integral  as 


V 


1  ■  ow*  2 l*  [d\  •••  4>r  r  £u  “• 


21) 


^(1+  1) 

where  y  indicates  summation  over  all  possible  single  cluster 

diagrams  which  can  be  found  from  (1+1)  given  particles. 


Many  of  these  cluster  integrals  can  be  further  simplified. 
Suppose  that  a  cluster  of  (£  +  1)  points  may  be  disconnected  by  breaking 
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it  at  a  single  point,  r  .  Then  b  may  be  written  as  the  product 

j  |+;i 

of  two  integrals,  one  corresponding  to  each  of  the  two  component  clusters; 
the  integration  is  carried  out  over  relative  coordinates  (r  -  r  )  ,  r 
being  considered  fixed;  in  principle  these  two  independent  integrals 
will  depend  on  the  fixed  position  of  r^,  relative  to  the  volume  boundaries. 
By  requiring  the  interaction  potential  to  be  finite,  and  by  having  V  go 
to  infinity,  this  dependence  is  eliminated. 

A  cluster  of  k  points,  in  addition  to  particle  1,  which 
cannot  be  reduced  by  a  single  break  is  termed  irreducible;  irreducible 
cluster  integrals  are  given  by 


^  r 

A"  £  ••• 

J 


3  3 

dr.  •••  d  r 
^  1. 


,1 


(4.22) 


^-i00 

where  y  indicates  summation  over  all  possible  irreducible  cluster 
diagrams  which  can  be  formed  with  (k  +  1)  given  particles.  For  k  =  0 
there  is  no  connection  between  point  1  and  any  other  point;  other  examples 


Irreducible  Integral  Diagram 


If  it  is  possible  to  define  other  even  more  specialised  classes  of  cluster 
integrals,  but  further  comments  about  this  will  be  deferred  until  the 
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discussion  of  the  cluster  expansion  for  the  pair  correlation  function 
*2 

C.  Cluster  Expansion  for  the  Pair  Correlation  Function 

In  this  section  a  cluster  expansion  for  the  pair  correlation 

function  and  for  the  potential  of  the  mean  force  will  be  considered. 

In  the  limiting  situation  where  both  the  particle  number  N  and  the  system 

volume  V  approach  infinity  while  the  particle  density  f*  -  H  remains 

constant,  g„  (r, ,  r_)  =  g„  (r,„).  From  Eq.  (4.16)  it  is  seen  that 
z  -z  Z  1Z 


-1  2-] 

82<r12)=QN  V 


1  f  r  .3  3 

J  3  3 


H 

r  S  °‘j 


(4.23) 


The  first  step  in  the  expansion  development  is  the  introduction 
of  the  "general  12  cluster”;  this  refers,  by  definition,  to  any  diagram 
in  which  a  number  of,  ,  particles  are  at  least  singly  connected  to 
either  or  both  of  the  fixed  points  1  and  2.  See,  for  example.  Fig.  6. 


/ 


Figure  6  Two  Examples  of  a  General  12  Cluster 


The  corresponding  cluster  integrals  are 


b.  <r12>  *  1  +  f12  *  ‘ 


b*  ‘V  -  jf?  •••  "V  2  Z<8>  TT'u 


(4.24) 
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where  y  refers  to  summation  over  all  possible  cluster  diagrams 
involving  in  addition  to  points  1  and  2, (other  points  which  will  now 
be  called  field-points.  The  contribution  to  the  integral  in  (4.7)  from 
a  diagram  containing  such  a  12-cluster  plus  other  disconnected  diagrams, 
factors  into  a  product  of  bj  (r^)  and  the  integrals  involving  the 
remaining  N  -  ((+2)  points. 

Many  of  the  general  12-cluster  integrals  can  be  reduced 
even  further.  We  define  a  "general  12 -irreducible  cluster"  with  k  given 
field-points  as  a  diagram  in  which  each  of  the  k  points  lies  on  at  least 
one  continuous  path  which  goes  from  point  1  to  point  2  without  passing 
any  point  more  than  once;  such  diagrams  for  k  =  1,  2,  and  3  are  given  in 
Fig.  7. 


(a)  <b)  (c) 


Figure  7  Examples  of  General  12-Irreduclble  Clusters 


An  even  more  restricted  classification  may  be  effected  if  the 
k  field-points  fall  into  two  or  more  groups  such  that  points  in  two  dif¬ 
ferent  groups  are  Interconnected  only  by  paths  which  pass  either  through 
point  1  or  2.  If  this  is  the  case,  the  diagram  splits  into  Independent 
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sub-diagrams  by  cutting  at  both  point  1  and  point  2.  The  corresponding 
general  12  cluster  integral  can  then  be  written  as  a  product  of  indepen¬ 
dent  integrals,  but  where  each  factor  is  a  function  of  r^>  A  diagram 
in  which  no  further  cutting  is  possible  is  called  "simple  12-lrreducible"; 
see,  for  example.  Fig.  8. 


k  =  1  k  =  2  k  =  3 
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A  -  £  J  \  *\  d\  f>3  £34  f42  [2+4  f14  +  £14  Q 

where  Che  numerical  factor  in  brackets  refer  to  the  number  of  diagrams 
in  which  only  the  ordering  of  field-particles  is  different. 

It  is  now  possible  to  write  down  an  expression  for  82^12^* 
it  requires  only  that  one  be  able  to  find  an  expansion  in  powers  of 
density,  ,  which  has  simple  1  2-irreducible  cluster  integrals  as 
coefficients.  Since  further  detail  is  concerned  only  with  counting  and 
collecting  the  various  contributions,  the  result  will  be  written  down 
immediately, ^ 

oO 

w2  (ri2)  =  “<r12)  *  S  P*  (r12V°  (4*27) 

"  *2  (ri2) 

*2  <r12)  *  *  (4.28) 

These  last  two  expressions  are  the  desired  result;  their 
application  to  the  case  of  the  long  range  potentials  Is  the  subject  of 
the  next  section. 

D.  The  "Debye-Chaln"  Type  of  Expansion ^ 

The  Irreducible  cluster  expansion  In  powers  of  the  particle 
density, ^  ,  for  the  pair  correlation  function,  gg  (r^g)  *nd  the  cor¬ 
responding  potential  of  the  mean  force,  Wg  (r^g),  was  the  subject  of  the 
previous  section.  It  was  Indicated  that  it  is  possible  to  label  each 
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irreducible  cluster  by  the  number  of  field-points  k,  end  the  number  of 
connecting  lines,  ^  ,  occurring  in  the  given  cluster.  For  e  specific 
k  value,  the  smallest  possible  value  of  is  given  by  k  +  1,  and  the 


largest 


by  [2  (k 


+  2)(k  +  1)  -  1 


ft. 


In  the  usual  expression  for  the 


irreducible  cluster  expansions,  as  shown  in  Eqs.  (4.27)  and  (4.28), 

the  diagrams  are  classified  solely  by  k,  all  values  of'ft  corresponding 

to  this  k  being  grouped  together.  When  the  particle  system  under 

consideration  involves  very  short  range  interactions,  i.e.,  neutral 

molecules,  there  is  no  trouble;  in  fact,  this  procedure  proves  quite 

satisfactory;  when  the  range  of  the  interaction  r  is  much  less  than 

-1/3  M 

the  average  particle  separation,  r  7C  P  •  Thus,  the  condition  for 

0  3 

the  "k  only"  grouping  to  be  valid  is  that  yOr^  «1.  Now  consider  such 

a  short  range  potential,  u^j,  which  is  large  and  repulsive  at  very  small 

distances,  slightly  attractive  at  the  potential  minimum,  and  drops  to 

sero  rapidly  for  r  J  r  .  In  this  instance  f  is  or  order  unity 

inside  the  region  r  ,  and  drops  to  zero  rapidly  outside;  the  con- 
R 

tributions  to  the  cluster  integrals  f3 ^  from  an  irreducible  diagram 

with  k  field  points  and  JJ  "lines”  is  then  of  order  ( p>  r^)^  ,  regardless 

of  the  JL  values;  here  we  are  excluding  factors  outside  the  integral, 

1 

such  as  vi  ,  etc. 

3 

An  effective  expansion  parameter  here  is  p  r^  ,  and  there  is 
no  need  to  regroup  any  terms  according  to  X  values.  In  the  event  that 
the  potential  u^  has  a  large  (attractive)  potential  minimum,  f  , 
becomes  large  and  positive  in  the  vicinity  of  this  minimum.  Thus,  for 
a  given  value  of  k  the  largest  contribution  to  jS ^  comes  from  diagrams 
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with  fairly  large  values.  Again,  there  is  no  reason  to  regroup  the 

expansion  terms.  The  situation  is  much  different  when  the  interaction 

potential  is  long  range;  here  the  "k  only"  arrangement  is  not  correct. 

In  accordance  with  both  classical  mechanics  and  chapter  three, 

a  plasma  model  will  now  be  adopted  which  involves  only  long-range, 

repulsive  interactions;  thus,  consideration  of  attractive  potentials  in 

general,  and  infinitely  attractive  potentials  in  particular  is  avoided. 

This  plasma  is  composed  of  singly  charged  ions  and  electrons,  and  involves 

interaction  potentials  (at  the  average  distance  of  particle  separation) 

-1/3 

that  are  small  compared  to  thermal  energy:  i.e.,  ui j (ri j  £  P  )  «  kT; 
this  condition  implies  a  high-temperature,  low-density  system. 

In  this  instance  one  hopes  to  find  another  small  parameter 
in  terms  of  which  the  cluster  expansion  may  be  constructed.  The  potential 
model  will  now  be  specified  in  three  distinct  regions: 

(1)  r  less  than  some  small  distance,  r ; 

H 

u^(r)  is  large  and  positive. 

(2)  rM  - r  -  x  * wher*  x  y?  »  utj 

is  small,  positive,  and  of  order  of  magnitude  £  , 
where  &  is  a  constant  much  less  than  1. 

*■3 

(3)  r  yy  X  ;  u^(r)  decreases  more  rapidly  than  r  . 

-u(r) 

It  is  further  assumed  that  (r^  //,)«  6  «  l.  The  function  f(r)  =  e  -1 
is  then  approximately  equal  to  -1  in  the  region  (1),  of  order  £  in  region 
(2),  and  is  negligible  in  (3).  Therefore,  the  total  volume  integral 
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(over  all  space)  is  of  order  6  X  Now  only  cases  where  Che  range  of 

Che  poCenCial  is  long  compared  Co  Che  average  parCicle  separation,  i.e. , 
*1/3  3 

^  y  yO  ,  or  X  ^  1  will  be  considered.  ReCuming  Co  Che  order 

of  magniCude  of  Che  irreducible  clusCer  integral,  ^  ^  ,  it  is  assumed 

Chat  each  integral  contains  k  field-points  with  coordinates  r  and 

5  V  3  k. 

"lines11.  The  order  of  magnitude  here  is  given  by  £  ( )  ,  or 

Q  o  k 

£  )  •  Although  £  is  by  definition  small,  Che  quantity 

^6-X"*  may  be  large  for  sufficiently  large  X  •  It  is  therefore  useful 
to  regroup  the  cluster  expansion  terms  for  g^  (r^)  according  to  the 
value  of  H  -  k;  the  summation  is  over  k  values  for  fixed  JL  -  k.  This 
procedure  is  known  as  the  cluster  expansion  in  complexities. 

Terms  in  Eq.  (4.27)  which  correspond  to  values  of  Jt  -k  =  1, 
correspond  to  the  "simple  chain"  diagrams  illustrated  below  in  Fig.  9. 
Recall  that  the  term  f  ^ ,  corresponding  to  k  =  0,  is  missing  from 
Eq.  (4.27). 


Figure  9  Clusters  of  the  "Simple  Debye  Chain"  Type 

The  sum,  8^  (r^),  of  the  terms  in  Eq.  (4.27)  with  Jt  -k  =  1  (k  £  1) 
is  given  by 

S,<r12>=2_,  PU  \  •**  W3'  •••  ***  •••  f(r  )  (4.29) 

1  n«l  /  J  J  3  n+2  13  tr+2,  2 
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Setting 


(4.30) 


it  is  observed  that  (r^) 
the  function  j  (r^)  ’  thusa 


is  the  Neumann-Liouville  expansion  for 


r«12>  -  *  **>  rJ4 « <ri3>  r*„> 

Taking  Fourier  transofrms,  we  find  that  the  resulting  equation  may  be 
solved  by  convolution  techniques,  yielding 

/“'(k)  =  f  (k)  -  fi(V.)  r 00 

or,  (4.32) 

r<k>  -  _ 

1  -  ft  (k) 


The  solution  of  Eq.  (4.31  or  Eq.  (4.32)  represents  in  closed  form,  the 

summation  of  all  diagrams  with  Jl  -  k  =  1,  for  avbitrary  k;  convergence 

(2) 

is  now  a  question.  It  can  be  shown'  '  that  in  the  range  - oo^ pf (k) <  +  1, 
Eq.  (4.32)  is  justified. 

If  only  the  X  =  k  =  1  terms  are  used  in  the  expansion,  we 
get  as  an  approximation  to  the  potential  of  the  mean  force,  w^: 

w2  (r12>  =  u  <rl2>  ”  Si  (r12> 


=  -r ;ri2)  +  u<rl2> +  f(ri2) 


(4.33) 
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If  u  *nd  f  are  of  order  £  «  1  throughout  the  range  of  the  potential 

.  2 
then  the  error  made  in  omitting  terms  with  Ji  -  k  =  2  is  of  order  €.  . 

Also  note  that 


=  u 

u  +  f  =  e  -  1  +  u 


t-f 


+ 


(4.34) 


2 

and  hence  is,  in  the  first  approximation,  of  order  £  .  Using  this 

first  approximation  f  is  replaced  by,  -u  ,  in  Eq.  (4.31)  giving 


-2°  <t12>  '  “  (t12> 


f 


j3 
d  r„ 


(rl3> 


(r32) 


(4.35) 


Applying  the  transform,  convolution  technique  to  this  equation, 


w  °  (k)  =  -  ~u  W - 

2  1  •  flu  (k) 


i-r 


If 


c  -cc  r 

-  (t„>  -  S-  .  12 

12  r12 


(4.36) 


s  =  ~ 

C  XkT 


where 

units  of  the  Debye  length,  X  .  Then, 


and  r^  *-s  measured  in 


-  u  (k)  = 

k 


/  i  I  =9<rl2 
dr i2  rl2  »in!kr12] *  r12  e 


4T€rA? 

k2  +<*2 


.37) 
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.(52) 


Now  it  can  be  shown ^  '  that  47T€^>  =  1.  Thus,  we  have 


w  (k)  = 


A.U.£L 


k2  +«2  +  l 


(4.38) 


oO 


W2°  (r12)  =  ‘L^~  \  dk  k  *in  (kr)  \-s - -1 

2  12  -rrr12  \  v  'Lk2+oc2  +  iJ 


ii  ri  , 

irr12  L  2  J 


(4.39) 


For  c<  =  o 


o  .  ,  6  “  ri2 

w  (r)  =  —  e  Li 

2  r12 


(4.40) 


which  is  just  the  Debye-Huckel  result.  If -  1, 


»2°(r) 


”  f*'  r12 


(4.41) 


Note  that  if  u  (r)  is  large  and  positive  over  even  a  snail 

region  of  space,  Eq.  (4.35)  becomes  a  very  poor  approximation  to  Eqs. 

(4.31)  and  (4.32).  However,  terms  corresponding  to  1  °  k  2  2  which  are 

2 

neglected  in  equation  (4.33)  are  still  of  order  £  .  This  implies  that 

for  small  £  they  involve  integrations  to  which  this  small  region  of 
space  contributes  little. 
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E.  Higher-Order  Terms 

Thus  far,  only  the  £  -  k  =  1  summation  has  been  discussed. 

Terms  for  X  1  may  also  be  included  as  a  sum  of  a  finite  number  of 

integrals  involving  the  f  (r)  previously  defined  in  Eq.  (4.30). 

Choose  a  typical  irreducible  cluster  diagram  occurring  in  Eq. 
(4.27)  with  arbitrary  values  of  k  and  such  as  the  diagram  shown  in 
Fig.  10-a. 


(a)  (b) 

Figure  10  An  Example  of  a  Simple  12-Irreducible 

Cluster,  In  Terms  of  "Wiggly  Line"  Diagrams. 

In  this  diagram,  each  of  the  numbered  points  3  to  7  are  connected  directly 
to  at  least  three  nearest  neighbor  points,  whereas  each  of  the  unnumbered 
points  are  connected  to  only  two  neighbors.  Fig.  10-b  again  shows  this 
cluster  diagram,  but  this  time  in  terms  of  "wiggly  lines",  where  each 
"wiggly"  line  represents  a  "simple  chain"  (Fig.  9).  How  all  diagrams 
which  correspond  to  the  same  "wiggly  line"  diagram  are  considered  together; 
although,  the  number  of  lines  or  links  making  up  each  "wiggly  line"  may 
be  different  in  each  case,  the  numerical  value  of  Jj  -  k  is  the  same  for 
all  these  diagrams. 
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In  way  of  summary,  the  following  "wiggly  line"  characteristics 
are  listed: 

1.  Each  field  point  lies  on  at  least  one  continuous  path 
between  points  1  and  2. 

2.  Any  pair  of  field  points  is  connected  by  at  least  one 
continuous  path  which  passes  neither  through  point  1  nor  2. 

3.  Each  field  point  is  directly  connected  to  at  least 
three  nearest  neighbor  points;  if  it  were  connected  to  only  two  nearest 
neighbors,  it  would  be  a  mere  link  in  the  chain  which  contradicts  the 
definition  of  a  "wiggly  line"  diagram. 

4.  Any  pair  of  points  may  be  joined  by  more  than  one 
wiggly  line. 

Hote  that  point  3.  is  in  general  not  applicable  to  the  usual  irreducible 
cluster  diagrams  and  that  point  4.  never  applies  to  such  diagrams. 

Finally,  an  examination  must  be  made  into  the  nature  of  the 
wiggly  line  factors  which  occur  in  the  integrands  of  the  summed  cluster 
integrals.  Consider  two  points,  (  i,  J  )  connected  by  a  single  wiggly 
line.  This  line,  »  represents  the  sum  of  the  direct 

link,  f  (r^)  ,  plus  all  the  simple  chains  with  an  arbitrary  number  of 
links  (see  Fig.  9).  But,  this  sum  is  singly  the  function  /”"*  sf+f 
defined  in  Eq.  (4.30).  The  general  rule  is  that  a  single  wiggly  line 
between  points  i  and  j  stand  for  a  factor  in  the  Integrand;  this  may 
be  compared  with  f  (r^)  for  the  ordinary  cluster  integrals.  We  next 
consider  n  parallel  connections  between  i  and  j;  here  at  most  one  of  the 
n  connections  can  stand  for  the  direct  link,  f  (r^).  The  factor  in 
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the  integrand  is  then  the  sum  of  two  terms 5  the  first  term  which  does 
not  contain  the  direct  link,  f  (r  )  ,  consists  of  the  product  of  n 
simple-chain  factors ,  C,  s  and  to  account  for  the  possible  chain  permu¬ 
tations,  a  factor  -  is  included.  The  second  term  contains  the  direct 


f-link,  and  in  addition,  n  -  1  simple  chain  factors.  For  each  n-fold 
wiggly  line  connection  between  i  and  j ,  there  is  a  factor  in  the  integral 
of  the  form, 


17 


<H,>  -  i  S°  <«„>  +  ^  Hri})  S°  1  (r^j )  (4.42) 


n: 


Having  thus  obtained  the  function  /""*  or  ^  from  Eq.  (4.31) 
or  (4.32)  we  have  an  explicit  integral  for  each  wiggly  line  diagram. 

The  value  of  £  -  k  is  equal  to  the  number  of  wiggly  lines  minus  the  number 
of  field  points;  ji  -  k  =  1  implies  a  single  such  diagram  while  jL  -  k  =  2 
gives  rise  to  three  graphs  as  seen  in  Fig.  11. 


Figure  11  The  Wiggly  Line  Diagrams  For  JL  -  k  =  2. 


Returning  to  the  potential  of  the  mean  force,  and  considering 
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terms  up  to  and  including  -  k  =  2,  we  find 


This  is  the  expression  we  have  been  seeking. 

In  order  to  make  some  comment  about  the  magnitude  of  successive 
terms  in  this  expansion,  we  again  consider  a  repulsive  potential  of 
range  X  ,  with  u  and  f  of  order  £ <<  1  (r  <  X),  but  with  |pcX3 /  >_  1. 

From  consideration  of  the  Fourier  transform  it  may  be  shown  that  j)  |  f(k)  | 
is  large  for  small  k  which  in  turn  implies  that  /"I '(k)  X  1  for  small  k. 

Thus,  it  might  be  expected  that  ^""*(r)  is  comparable  with  f  (r)  •*-  €  for 
small  enough  r,  but  decreases  more  rapidly  than  f  (r)  for  r  yjX  . 

|r<k)|  for  small  k  is  of  order  ^  d\  (r)  «-  £  X  ^  and  we  therefore 
further  expect  that  (^€X  )  will  be  of  order  unity.  If  this  is  the  case, 
the  order  of  magnitude  of  any  term  for  any  value  of  jl  -  k  may  be  estimated 
to  be 


4  -  k  .  k  JL  -  k 


(4.44) 


Hence,  terms  corresponding  to  large  complexity  —  (1  -  k)  =  large  •—  may 
be  expected  to  contribute  little  to  the  cluster  expansions. 


V. 


CORRECTION  TO  THE  DEBYE -HUCKEL  THEORY  AND  ITS  APPLICATION 
A.  Introduction 

In  this  chapter  the  theory  developed  in  the  preceding  chapter 
will  be  used  to  derive  an  expression  for  the  potential  of  the  mean  force, 
Wg,  which  is  then  applicable  to  either  the  Coulomb  or  shielded  potential. 
Correction  to  the  usual  Debye-Huckel  expression  for  w^  will  be  given  to 
second  order  in  complexity.  Furthermore,  the  results  of  these  corrections 
will  be  scrutinized  in  light  of  the  discussion  in  Chapters  II  and  III. 

Corrections  of  the  nature  proposed  have  previously  been  calcu- 
(52) 

lated  by  D.  L.  Bowers  for  the  case  of  the  Coulomb  potential;  however, 
in  his  work  he  made  the  error  first  of  using  inaccurate  tables,  and 
secondly,  of  using  an  incorrect  computer  technique.  Since  his  calcula¬ 
tions  and  detailed  data  are  not  available  in  the  open  literautre,  it  is 
not  possible  to  determine  the  exact  importance  of  these  errors.  In  view 
of  this  uncertaintly  it  was  felt  worthwhile  to  repeat  the  Coulomb  calcu¬ 
lation  in  addition  to  considering  the  shielded  field  case. 

The  expression  for  the  potential  of  the  mean  force,  w^,  will 
be  evaluated  using  first  only  the  -  k  =  1  term  and  then  including  the 
X  -  k  =  2  contribution.  Thus,  we  must  consider 


where 


w2(ri2)  =  «(r12)  +  f(r12)  -  P (r^)  -  X(r12>,  (5*1) 


,  ^  _  U  <r12> 

u<ri2)  =  - 


kT 


“u(ri J 

f (r12>  =  e  12  -1 


r (ri2>  =  f(r12>  +P 
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Mote  that  £  gives  a  measure  of  the  strength  of  the  interaction  energy 
relative  to  the  thermal  energy. 

As  was  shown  in  Chapter  IV,  in  the  event  £  approaches  zero,  H 


approaches  which  is  just  the  Debye=Huckel  result.  In  the  present 

o 

notation, 


(5.3) 


The  valuesof  €  which  will  be  considered  are  £  =  0.1  and  0.2.  These  are 

considered  small  values;  from  Eq.  (4.54),  it  is  known  that  the  order  of 

magnitude  of  the  ratio  of  the  A  -  k  —  2  t  erms  to  those  for  which  ^ -  k  =  1 

will  be  £.  ,  while  the  ratio  of  higher  order  terms  to  the  (  %  ~  k)  =  1 
2 

terms  will  be  £  •  Therefore,  it  is  expected  that  the  Jl  -  k  =  2  terms 

will  provide  a  satisfactory  stopping  point. 

To  simplify  notation,  we  will  set  r^  =  r  in  the  remainder 
of  this  chapter. 

B.  Evaluation  of  /~* Cr) 

From  the  theory  of  Fourier  integrals  we  may  write  the  fol¬ 
lowing  3-dimensional  transforms. 


aO 

P  (k)  -  42T  r  d 

k  J 

o 

oO 


sis(kr)F  (r) 


(5.4) 


F  (r)  = 


2  TT2t 


dk  k  sin  (kr)  F  (k) 


The  Fourier  transform  f  (k)  has  been  shown  in  Eq.  (4.32)  to  be  simply 
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“  u  (r) 

Since  f  (r)  =  e 
transform  of  f  (r)  as 

f  (k)  =  ilE 

k 


€ 

*  X.  r 

-  1  and  u  (r)  =  — 

e  we  may  write  the 

r 

oC 

r\ 

r 

c  -  X  r 

~  —I  p  -n 

dr  r  sin  (kr) 

e  r  -  l]  (5.5) 

It  is  possible  to  effect  an  approximate  solution  of  this  integral  by 

brute  force  techniques ,  but  the  solution  for  (X-^  0  is  so  cumbersome 

that  it  is  of  little  use.  In  the  event  that  one  is  interested  in  the 

limiting  case  of  £<  — *-  0  ,  the  resulting  solution  may  be  expressed  in 

terms  of  ker  functions: 

2 


f  (k)  -  -  .  2  ker2m 

«=0  k2 

where  (5.6) 

„  1/2 

=  (4  £  k) 


If,  however,  it  is  desired  to  solve  the  expression  for  f  (k), 
for  arbitrary  c<  ,  another  method  must  be  devised.  An  approach  which 
avoids  the  difficulties  mentioned  before  is  as  follows: 


r 


Let  ^  e  =  R  .  The,  the  expression  ^e 

can  be  expressed  in  terms  of  Mellin  transforms. 


•R 


-R 


e.  =1 


<T+  io& 


ds 


i  60 


-i<<r<o  (5.7) 
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Inserting  the  correct  form  for  R  and  substituting  in  Eq.  (5.5),  we  find 


(T+t00 


fOO- 


^  dr  r  sin  (kr)  /^(s)^  ^  e  J  ds 


c r-  i  “5 


Simplifying  and  inverting  the  order  of  integration,  we  have 


<rt/w>  o 

f  f 

f(k)=  *  2¥i  r<«>€ 

cr-ioo  o 


-s  s+1  ocsr 

s)£  ds  dr  r  sin  (kr)  e  (5.8) 


Now  define: 


t  *  r/g  r  *6  t 


*>£=a  ^  <xr  =  °<€  t  =  at 


Using  these  definitions,  we  get 


t8  e88t 


r  sin(kr)dr  =  g2  t  sin(k£  t)  dt 


Thus, 


(5.9) 


2 

=  €  t  sin  (bt)  dt  b  =  -  k 


C3“+  (oO  CC 

r  r 

s+1  ast 

•  — - —  /  (s)  ds  dt  t  sin(bt)e 

2  7T  i 

cr-i«o  o 


(5.10) 
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This  integral  may  be  solved  by  contour  integration.  We  choose 
the  contour  as  shown  in  Fig.  12. 


Figure  12 


68 


By  Cauchy's  residue  theorems, 


C  C 

1  2 


=  2  IT  i  y  Residues 


(5.13) 


It  is  shown  in  Appendix  A  that  the  contribution  from  =  0.  Thus, 


CT+- 1  00 


s)  [* (s+2)  sin  ["  (s+2)  tan  1 


2  +  b2  2 


=  2TiZ 


Residues  (5.14) 


cr-i«o 


Included  inside  the  chosen  contour  are  poles  at  all  integer  values  on  the 
negative  real  axis;  these  are  due  to  singularities  in  the  product  of  gamma 
functions.  At  s  -  -1  there  is  a  single  pole  while  for  s  £  -2,  the 
poles  are  all  double.  At  s  —  -1  the  residue  is 


Rgp  =  Residue,  single  pole 


,  T  -1  bl 
sin  L  tan  j  J 

[777p 


(-1)  P (®t2)  aLn  l  (s+2)  c.n  '( j/ 

o4-9  r  (5. 15) 


[->  “if 
(as)2  +  b2J  2 


£  G*2  +  kzl 


For  all  other  poles. 


R  =  Residue,  double 
DP 


P°le’  =  [  A  {  (S4*>2  ^ iS)  ^  (s+2)F(s)]]  (5*16) 
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where 


F<s)  * 


sin  (s+2)  tan  (~Vl 
[(,s)2  +  b2]  ^ 


i  —  4LC si) 

P (n-l)pin+l)  _  ds 


_  Fi*,n 


_  KmJ 

F  (n+l)  i 


^(n-l)  F  (n+l) 


(5.17) 


Performing  the  operation 


d  F  (b) 


yields, 


d  F  (s) 


L.,=K+b^ 

+  sin  X  In  ((an)2+  b2)  '  +  **  02*W"j*-  ? 

(an)  t  bz  J  J 


«  f  -  (n-2)  tan  ( -M 1 
»-  an  J 


where 


Using  this  expression  for  fLF_I®2  |  in  gq.  (5.17)  gives  for 

d  s  s  =  =n 

the  residue  at  a  double  pole. 


R  = 
DP 


[(an)2  +  b*P 

F (n+l)  Ft n-D 


)  Tab  (n-2)  -1  .  t  .1 

■( C0‘  *  L  -  -  (in 
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+  sin 


1/2  2 

in  X  T  In  ((an)2+  b2)  +  a— n  in"2) 

L  (an)2  + 


(_>*  +  **  <■*»]] 


(5.18) 


where 


'V  (n)  »  JZini 

r  (n) 


Substituting  Eq.  (5.15)  and  (5.18)  into  Eq.  (5.14)  one  finds  that 
(5.12)  becomes 

3 


f  (k)  =  4H3L 

b 


=  42V 
b 


I 


Residues 


y  [~(an)2+  b2] 

(.2  +  b2)  1/2  r(n+l)/]n-l) 


n-2 
2  ... 


(5.19) 


where 


n 

j 

) 

ii 

cos  X 

1 

LJ 

f  — b— (n_2)  _  tan‘1/JL')l 

(an)2+b2  an  J 


+ st°  *  [in((“)2 + b2>l/2 + -  r<««> 

(An)  +b  -J| 


If  0<  -  0,  and  if  the  definitions  in  Eq.  (5.9)  are  recalled,  this  general 
expression  reduces  to 


[eo 

-I 

n=2 


n-1 

f  (k)  — SJK.  )  1+  >  - jjz 

k4 


2  CMl  2 


i/(n+l)/(n-l)l 
+  sin  [  ln(e  k)-An+l)  -  f(n-l)jJ 


(5.20) 


As  is  shown  in  Appendix  B,  this  may  be  expressed  in  terms  c£  modified 
Bessel  functions: 
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f  (k)  =  -  •  2  ker 


*  I’ 


1/2 


=  (4€k) 


(5.21) 


Poro<  =  1,  corresponding  to  the  Debye-Huckel  shielded  field,  f  (k) 


becomes 


f  (k) 


k2  +  1 


oo  2  2 


1  — 


k  +1 


I 


C»V) 


n-2 

2 


n-2  p (n+1)  /~(n-l) 
=  (cos  [(n-2)  tan  (  J)]  [  tan  J)  - 


(5.22) 


+  sin  [(n-2)tan‘1(J)][ln;€(n2+k2)  *+  f  (n+l)- ^(„-l) 

n  *  k  _ 


It  is  Interesting  to  note  that  in  Eqs.  (5.20)  and  (5. 22), the 

first  term,  due  in  each  case,  to  the  single  order  pole,  corresponds  to 

the  Debye-Huckel  result:  if  only  this  term  is  kept  and  the  analysis 

carried  on  one  finds  that  /”*  =  '  =—• w  .  The  fact  that  all  other 

o  2 

terms  arise  from  double  poles  is  interpreted  in  the  following  manner: 

In  the  first  approximation  one  takes  only  the  first  term  in  the  expres¬ 
sion.  If  higher  order  corrections  are  desired,  it  is  necessary  to 
include  the  entire  summation  over  double  poles  and  not  just  the  first 
few  such  terms.  Actual  calculations  confirm  this  last  statement. 

The  above  analytic  expressions  allow  one  to  determine  /"^(k). 

However,  it  is  not  possible  to  get  the  inverse  transform,  /"*( r),  analytically. 
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is  solved  by  computer  techniques.  Pertinent  information  about  the 
computer  approach  is  contained  in  Appendix  D.  It  should  be  mentioned 
here,  however,  that  in  both  instances  for  which  numerical  work  was 
done,  o <  =0,  1,  the  needed  values  of  f  (k)  were  generated  by  the 
machine  from  Eqs.  (5.20)  and  (5.22)  and  are  not  tabular  values.  The 
above  integration  was  performed  for  £  a  0.1,  0.2  for  both  =  0  and 
o£  =  1;  the  results  are  listed  in  Table  I.  Once  /"*( r)  is  known, 
it  la  immediately  possible  to  find  expressions  for  w^(r)  and  g2(r): 


w2(r)  ■=  u  (r)  +  f  (r>  -  T  (r) 
-  w  (r) 

g2  (r)  ~  e 


The  expressions  found  for  w2(r)  by  this  method  are  compared  graphically 

with  the  corresponding  DebyeHuckel  results,  (r)  =  w  °(r)  ,  in 

o  ‘ 

Figs.  13,  14,  15,  and  16.  It  should  be  pointed  out  here  that  g  can 

n 

be  related  to  the  function  s^  discussed  in  Chapter  III^*®);  specifically, 
setting  n  *  2,  we  find  for  a  fluid  system  that 

g^(r)  =  1  +  s2(r) 

It  is  interesting  to  note  that  the  corrected  curve  for 
w2  £  6  =  0. 1  j  «*  =  1 ?  lies  just  above  the  Debye-Huckel  curve,  while 
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result  for  f  €  =■  O.TLj®!  =  ljli.es  considerably  below 


the  cor- 


responding  Debye-=Huckel  curve.  This  implies  that  in  the  first  instance 
the  Debye-Huckel  approximation  overestimates  the  shielding,  while  in 
the  second  case,  it  underestimates  the  shielding.  That  such  a  result 
is  at  least  possible  may  be  seen  by  considering  the  graphs  of  f  (k) 
for  various  values  of  £  as  shown  in  Fig.  17.  The  Debye-Huckel 
result  corresponds  to  £  =  0  ,  which  fact  implies  that  it  is  valid 
only  in  the  infinite  temperature  limit,  T  — ►  .  Now  observe  that 

as  the  magnitude  of  £  increases,  the  slope  of  the  curve  for  smaller 
values  of  k  increases  noticeably,  while  the  positive  region  of  the  curve 
becomes  much  larger  and  shifts  toward  the  smaller  k  values.  From  actual 
computer  runs,  it  is  known  that  k  values  above  150  make  negligible 
contributions  to  the  integral  expression  for  f  (r)  .  As  £  is  increased, 
the  curves  for  f  (k)  change  as  indicated  above,  but  the  positive  portion 
of  the  curve  does  not  become  important  until  €  v'O.l;  the  change  in  the 
curve  which  is  important  below  £ ~  0. I  is  the  alteration  of  slope.  This 
importance  is  apparent  when  one  realizes  that  this  function  is  multiplied 
by  sin  (krj  in  the  Integrand.  Thus,  the  slope  change  in  the  region 
below  £  =0.1  is  the  primary  cause  of  deviation  from  the  Debye-Huckel 
curve.  When  the  positive  segment  of  the  f  (k)  curve  becomes  appreciable, 
in  the  range  c £  significant  k  values,  this  feature  of  the  curve  becomes 


the  more  important  factor,  and  the  corrected  expression  for  w^ir)  lies 
below  its  Debye-Huefcel  counterpart.  This  brief  discussion  is  designed 
only  to  indicate  that  such  an  £  dependent  shift  is  plausible. 
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TABLE  I 
r (r)  vs.  r 


0 


£  =  o.i 

r"(r) 


0.0 

=1.33390 

0.1 

-0.83742 

0.5 

-0.10685 

1.0 

-0.03620 

1.5 

-0.01523 

2.0 

-0.00705 

3.0 

-0.00177 

€  -  o.i 

T(r) 


0.0 

-0.99404 

0.1 

-0.51042 

0.5 

-0.09401 

1.0 

-0.02454 

1.5 

-0.00853 

2.0 

-0.00304 

3.0 

-0.00049 

€  =  0.2 


r 

rw 

0.0 

-0.82915 

0.1 

-0. 68575 

0.5 

-0.16869 

1.0 

-0.04958 

1.5 

-0.01857 

2.0 

-0.00681 

3.0 

-0.00125 

6  =  0.2 


r 

rw 

0.0 

-0.91017 

0.1 

-0.62486 

0.5 

-0. 15745 

1.0 

-0.03720 

1.5 

-0.01233 

2.0 

-0.00377 

3.0 

-0.00040 

Constants  employed  in  these  calculations  were: 
IT  -  3.14159265 

lj  =  Euler's  Constant  =  .577215665 


Figure  13 

w2  (•*)  vs  r  for 

a  =  1.0,  «  =0.1 

I  -  k  =  1 


Debye-Hiickel 


Corrected 
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000 


Region  of  Crossover 
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For  the  cases  OC=0,£  =0.1,  0.2,  the  corrected  w^(r) 

#* 

curves  lie  below  the  corresponding  Debye-Huckel  graphs  until  r  approaches 
unity;  in  this  region  the  corrected  curves  cross  the  approximate  curves 
and  remain  above  as  r  is  increased.  The  w^(r)  curve»  for  €  -  0.2, 
deviates  more  from  the  Debye-Huekel  result  than  does  the  £—0.1 
figure;  furthermore,  the  crossing  point  in  this  case  occurs  for  larger 
r  values  than  is  the  case  with  £  =0.1.  See  Figs.  15  and  16.  This 
implies  that  for  r  <  1  the  Debye-Huckel  approximation  overestimates 
the  shielding,  while  for  r  ^  1  it  underestimates  it. 

Now  consider  the  contribution  of  the  £  -  k  =  2  term  which 
are  contained  in  X  (r).  From  the  preceding  chapter,  we  use  Eq.  (4.  1). 


Also,  from  arguments  presented  In  the  previous  chapter,  it  Is  expected 

2 

that  the  magnitude  of  this  term  will  be  of  order  £  .  Now,  since  the 

expressions  for  /”* ( r)  derived  for  -  k  -  1  do  not  vary  greatly  from  the 

Debye-Huckel  results,  X  (r)  will  be  evaluated  setting  /”* (r)s  j  (r)  ; 

o 

this  implies  that  X  (r)  is  calculated  to  first  order  in  £  .  With  this 
substitution,  X  (r)  may  be  evaluated  analytically;  this  calculation 
for  arbitrary  oC  is  shown  in  detail  in  Appendix  C.  The  final  expression 
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ForCs^^  1,  the  expression  (5.24)  becomes 


w-sz  l;'17* 
16  ' 


/  -(2-/2)r  \  L  /  -/Tr 

2-<r  *  Tf?  e  _1 


(5.26) 


<r  +  ^)  (  E_(3/Fr)  -  E_  ((24^)r)  ) 


!  -<2+f2’)r  -3/T; 

*  ">  o  •  _ • _  — 


2+vt 


3f r 


These  functions  are  shown  graphically  for  both  values  of  oC  in  Pigs. 

18,  19  and  20,  21.  Note  that  the  contribution  of  the  J£_  -  k  =  2  terms 
is  much  less  than  that  from  the  corresponding  R  =  k  =  1  expression;  this 
is  especially  true  in  the  shielded  field  case. 

In  the  next  chapter  further  implications  of  the  results  dis¬ 
closed  in  the  present  chapter  will  be  discussed. 
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VI.  CONCLUSIONS 


Chapters  II,  III,  IV,  have  been  designed  to  show  the  interrelation 
between  the  spectral  line  profile,  the  electric  microfield  distribution, 
and  the  pair  correlation  functions  or  potentials  of  the  mean  force;  it 
was  seen  how  the  correlation  function  influences  the  microfield  distri¬ 
bution  which  in  turn  modifies  the  spectral  line  shape. 

In  Chapter  V  the  method  of  correcting  the  Debye-Huckel  pair  cor¬ 
relation  function  was  indicated  and  some  aspects  of  the  results  were 
disclosed. 

Ue  want  to  work  in  reverse  in  the  present  chapter,  that  is,  to  con¬ 
sider  first  the  implications  of  Chapter  V  and  then  to  apply  these  to  a 
discussion  of  the  preceding  chapters.  For  this  purpose,  it  is  found 
most  convenient  to  consider  tvo-body  effects  expressed  in  terms  of  the 
potential  of  the  mean  force  rather  than  the  pair  correlation  function; 
this  is  due  to  the  fact  that  in  Chapter  III,  it  is  this  quantity,  con¬ 
tained  explicitly  in  ^  ^  (y),  which  finally  enters  the  microfield 
calculation.  First,  the  potential  of  the  mean  force  will  be  considered 
for  the  Coulomb  example;  then,  due  to  its  direct  application  to  the 
quasi-static  aspects  of  the  line  broadening  theory,  the  case  of  the 
shielded  potential  will  be  discussed  in  some  detail. 

For  the  Coulomb  case,  the  corrected  expression  for  the  potential 
of  the  mean  force  indicates  that  the  Debye-Huckel  result  overestimates 
the  effective  repulsive  potential  for  smaller  values  of  r,  while  under¬ 
estimating  it  at  the  larger  r  values:  Close  interactions  are  more 
favored,  weak  ones,  less  so.  From  Figs.  20  and  21,  it  is  seen  that  the 
contribution  to  WgCr)  over  the  entire  r  range  from  the  JL  -  k  =  2  term 
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is  much  smaller  than  that  due  to  the  Jl  -  k  =  1  term.  Even  though 
rigorous  proof  is  missing,  the  supposition  that  higher  complexity 
terms  (L  -  k  2.  3)  will  be  negligible  seems  reasonable. 

The  £  =0.2  curve  differs  from  that  for  £  =  0.1  in  two  respects: 

(1)  the  magnitude  of  its  deviation  from  the  Debye-Huckel  graph  is  larger 
over  most  of  the  range  of  r  values,  the  region  of  crossover  being  the 
only  exception;  (2)  the  point  at  which  it  crosses  the  Debye-Huckel 
curve  occurs  at  a  larger  value  of  r. 

The  shielded-field  case  provides  results  which  are  more  interesting 
than  those  found  in  the  Coulomb  case.  A  comparison  of  Fig.  13  and  14 
shows  that  the  deviations  of  the  £  =  0.1,  0.2  cases  from  the  cor¬ 
responding  Debye-Huckel  results  are  not  only  different  in  magnitude, 
but  also  in  direction;  as  the  parameter  £  is  varied  from  0.1  to  0.2, 
the  relative  position  of  the  Debye-Huckel  graph  and  the  corrected  curve 
reverse.  For  £  =  . 1  the  corrected  curve  lies  above  the  Debye-Huckel 
figure  thus  implying  that  the  latter  result  underestimates  the  shielding 
effects.  When  £  =  .2,  the  correction  lies  below  the  Debye-Huckel 
curve,  implying  an  overestlmation  of  the  shielding.  For  neither  value 
of  £  does  the  correction  lie  equidistant  from  the  Debye-Huckel 
expression:  the  percent  difference  between  the  graphs  is  a  function 
of  r. 

Now  let  us  consider  how  this  behavior  of  the  shielded  field  might 
affect  the  ionic  microfield  distribution.  As  was  indicated  in  Chapter  111, 
the  low-frequency  microfield  distribution  is  strongly  dependent  upon  the 
nature  of  the  pair  correlation  function  or  the  potential  of  the  mean  force, 
even  though  these  quantities  are  much  smeller  than  their  companion  one-body 
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terms.  Mote  (see  Fig.  3)  that  the  effect  of  including  Debye-Huekel  pair 
correlation  is  to  shift  the  most  probable  field  strength  markedly  toward 
weaker  fields ,  while  simultaneously  increasing  the  magnitude  of  this  peak. 

It  is  now  possible  to  conjecture  about  the  effect  of  the  corrected 
potentials  of  the  mean  force  on  the  microfield  distribution.  In  the 
€  =  .1  case,  the  new  potential  is  everywhere  greater,  which  means 
that  all  possible  fields  are  slightly  less  probable  than  indicated  by 
the  Debye-Huckel  approximation.  Therefore,  it  would  be  expected  that 
the  corrected  microfield  distribution  would  be  shifted  little  if  any, 

For  £  =  .2,  the  relative  correction  becomes  larger,  as  r  is  increased; 
thus,  although  larger  fields  become  a  little  more  likely,  the  weaker  fields 
are  much  more  probable.  The  microfield  distribution  would  have  its  peak 
slightly  lowered  and  shifted  noticeably  toward  weaker  fields,  while  the 
strong  field  wing  would  be  raised  a  small  amount. 

It  is  fully  recognized  that  the  above  discussion  is  only  supposition, 
and  that  the  microfield  calculation  must  be  carried  out  for  various 
values  before  any  definite  statements  can  be  made.  However,  the  need 
for  such  a  revised  calculation  has  been  demonstrated. 

It  Is  interesting  to  consider  further  whether  it  is  possible  to  verify 
or  disprove  a  given  set  of  microfield  calculations.  One  such  possibility 
exists  in  the  case  of  the  low-frequency  or  ionic  microfield  distribution; 
this  involves  the  shift  of  spectral  line  profiles  as  shown  in  Chapter  11. 

In  order  tc  separate  the  microfield  effect  from  the  ocher  shift -inducing 
possibilities,  we  must  choose  a  spectral  line  for  consideration  where  the 
desired  effect  would  be  expected  to  be  large.  With  this  requirement  in 
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mind,  it  was  first  proposed  that  the  shift  of  the  peaks  of  the  Lyman  -jS 
line  be  examined  experimentally,  and  the  results  compared  with  those 
predicted  by  theory;  a  typical  Lyman  Q  profile  is  shown  in  Fig.  22. 


J  <U>) 


Figure  22  Typical  Lyman  j3  Profile 

The  experimental  study  of  the  Lyman  line  has  recently  been  done 

(45) 

for  a  given  temperature  and  density  condition  by  Elton.  However,  his 

results  are  not  conclusive  since  the  position  of  the  line  peaks,  strongly 
affected  by  self-absorption,  do  not  allow  one  to  discriminate  between 
currently  proposed  distributions.  There  is  another  possibility  which 
also  concerns  the  Lyman  @  line;  in  this  case  the  line  profile  would  be 
determined  holding  T  relatively  constant,  for  several  values  of  density. 
In  this  fashion,  one  should  be  able  to  observe  the  trend  in  the  peak  shift. 
Such  an  experiment  is  definitely  possible. 

It  has  been  established  that  the  correction  of  the  pair  correlation 
function  and  corresponding  potential  of  the  mean  force  should  produce  a 
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considerable  change  in  the  electric  microfield  distribution  function  cal¬ 
culated  using  the  Debye-Huckel  expressions.  This  change  will  result  in 
the  alteration  of  spectral  line  shifts,  which,  in  certain  cases,  are  subject 
to  experimental  verification. 


APPENDIX  A 


In  Che  evaluation  of  Che  integral  in  Eq. 


Chat  Che  contribution  of  \  was  zero,  i. 

J  c2 


(5.13),  it  was  stated 
e. , 
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which  for  large  Re  s 


=  sin 


[-  JL  +  i 
L  a  3 


bJ  +  ,,, 

s2 


For  Re  s 


sin  £  (s+2)  tan"^(--^— )J 


-  sin  _L 

a 


From  Sterling's  Theorem, 


ln/”(s)  =  ln(s-l)!  =  (s+-L)  ln(s-l)  -  (s-l)  +  lnflTp 


In  f( s+2)  =  ln(s+l)!  =  (s+3/2)  ln(s+l/2)  -  (s+1)  +  ln^TP 


Therefore, 


In  [/Is)  f  s+2)]  = 


(s-i-)  ln(s-l)  +  (s-hi)  In  (s+1) 


-  2s  +  In  2 r 


=  (2s+l)  In  s  -  2s  +  In  2TT 


(A.4) 


P(»)  As+2>  =2  TT  s2s  +  1  e‘2s 


Also,  for  large  Re  s 


[(«)2  +  b2  ]  ^  -  (,s) 


(A.  6) 
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This  implies  that  for  large  s, 


r<s)  r 'MD 

s+2 

[(as)2+b2j  2 


27T  s2s+2  e~2s 

s+2  s+2 

a  s° 


2  7Ts8  e-28 

-s  +  2 


(A.  7) 


Including  the  sine  term  we  find  that  the  product  of  s  times  the 
integrand  for  large  values  of  s  becomes 


s  *  Integrand 


=  2  7T  a8  e~2s  sin  (-  -gQ 

s  +  2 


(A.  8) 


This  is  in  general  not  equal  to  zero,  which  is  what  must  be  shown. 

After  a  change  of  variables  from  sto  |  =  s  -  |<rj  we 
determine  the  log  of  s  •  I  and  take  the  real  part.  This  must  be  equal 

to  —  c*C>  if  lim  V  - ^  0  . 

06  c2 


In  [  J  =  21T  +  (^+( T)  In  (^+<T)  -  2  (^+<T) 

-  q  +  CT+  2)  In  a  -  In  sin  (-  -L) 


=  In  21T  +  q  +<T)  In  ^  -  2^ 

-  (O'  +  2  )  In  a  -  $  In  a 


-  20“ 
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Let 

^  =  R  e*^  =  R  £  cos  <p  +  1  sin  <f>\ 

Then, 

RS  ln  L  J  *  [ln  2T1”  "  2<r  "  ^<T*+  2>  ln  *J 

+  (R  cos  ^  +  (T)  ln  R  -  R  <p  sin  <p 
-  2  R  cos  Cf  -  R  cos  ^  ln  a 


For  R  very  large,  this  reduces  to 


Re  ln  jjs  •  ij  =  R  ln  R  cos  (p  -  R  <f>  a±n(jO 


Subject  to  the  restriction  JIL  ^  (p 


lim 
R  - 


Re  ln 


[-*] 


oo 


(A.9  ) 


(A. 10) 


But,  this  is  just  the  region  of  integration;  hence,  Eq.  (A.l)  is  shown 
to  be  valid. 
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APPENDIX  B 

In  Chapter  V  it  was  stated  that  Eq.  (5.19)  reduced  to  Bowers9 
expression  for  Oi.  =  0.  This  will  be  demonstrated  in  this  appendix. 
Setting  o*.  =  0  in  Eq.  (5.19)  yields, 


f  (k) 


[oO 

i  +  Y . (fr.KI— _ 

n=2  r C“+l)  /W> 


(5.20  ) 


where 


(Sf)+  sin(Bf)  [lnek  -  ^(n+1)  -  f  (n-1)] 


In  order  to  prove  the  above  assertion.  Bowers'  result  will  be  shown  to 
be  identical  to  Eq*.  (5.20).  From  Chapter  V,  Bowers'  result  is 


f  (k)  = 


47Te 

,.2 


ker2  (^)  ;*£=(4€k) 


(5.21) 


By  definition, 


ker2  =  “  Re 


(B.l) 
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where  K  is  a  modified  Bessel  function  which  is  defined  as 
2 


oO 
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where 


JT 

2 


cos(lX)+  sin^JQZ^jln  (E  k)  -^(r+1)  -  ^(r-l)| 


which  is  identical  to  the  result  derived  from  Eq.  (5.19)  by  setting 0<  =  0. 


APPENDIX  C 


Solution  of  the  expression 


r 


<'l2>  'S’ 


dr. 


\T ^r13^ 


2  2 
r  <r32>  "  f  <r32> 


+  EL 
2 


dr3  dr4/(r13l  f<r34)  -  f2(r34^/  (r42> 


Jr <’ 


(4.43) 


using  the  Debye-Huckel  expression  for  A**  ,  and  letting  f  =  f  , 

o 


Ao  ~Tm  -  - 

o  r 

f  (r)  =  f  (r)  «  -  e"*r  (5*3) 

o  r 


It  is  here  convenient  to  introduce  bipolar  coordinates,  see  Fig. (C.l). 

y 


Figure  (C.l)  Bipolar  Coordinates 
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3 

d 


r3=  2TT  y  dxdy 


(C.l) 


Transforming  to  j,  r  coordinates,  we  use  the  relation 


x  X 

s  s 

dsdr  =  dxdy  (C.2) 

X-R  £ 

r  r 


Therefore, 

d3r  =  21T  y  dxdy  =  Ml*,  da  dr  (C.3) 

3  R 


First,  consider  the  evaluation  of 


Before  proceeding  any  further,  it  should  be  noted  that  as  in  the 
text,  the  interparticle  distances  are  given  in  terns  of  the  Debye  length,  X 


101 


A  consequence  of  this  Is  that  for  every  volume  element  appearing  In  the 
Integrals  to  be  evaluated,  a  factor  A  3  must  be  Included.  Furthermore, 
as  a  result  of  the  definitions  (5.2)  in  the  text, 


e  =  - — i - 

47T^A3 

or  (C.4) 

47T/°£A3  =  1 


Using  the  preceding  information  and  definitions,  one  has 


!  =  .,..(47T£g.A 

2  R 


oO 


(C.5) 
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R 


R 


If  one  pursues  this  calculation  further  and  in  doing  so  introduces  the 
functions 

O O 

(R)  =  f  - -  dr  (C. 7) 


R 

nO 
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Now  the  second  integral  to  be  evaluated  is,  II, 


Pl 

II  =  L-— 

2 


n 

A*r<v  dSr<'n)[r 


2  2 
’<*,,>  -  f  (r. 


It  will  be  noticed  that  the  second  integral  in  this  expression  is  equal 
to  I  (r^)  I  >  80 


r 

11  I  (r14>r >  <r42)  d\ 
_ 


Again,  using  bipolar  coordinates. 


(C.10) 


d  r4  =  2.Tir  d8  dr 


(C.ll) 


Therefore, 


-  I  (s)  sds  p  (r)  r  dr 


(C. 12) 


Making  the  appropriate  substitutions  and  carrying  out  the  rather  complicated 
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calculation,  an  expression  for  II  (R)  is  determined: 
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APPENDIX  B 


The  machine  calculations  for  the  work  reported  here  were  performed 
on  an  LGP-30  computer. 

The  numerical  work  may  be  separated  into  two  steps:  (1)  the  calcu¬ 
lation  of  f  (k),  and  (2)  the  evaluation  of  /"""( r). 

1.  Calculation  of  f  (k) 

From  equation  Chapter  V, 


f  (k)  - 


4  ire 


r  ** 

-X 


n-1 


IOI 


r  (n+l)/l [n-l) 


] 


(5.20) 


>  = 


IT 


cos(BX')+  sin(SX)|i!r1ie  k)  =  ^(n+1)-  ^(n-1)] 


This  is  the  expression  derived  from  equation  (5.. 19)  for c<  -  0;  the  remarks 
to  be  made  in  this  appendix  will  also  apply  to  the  more  general  result. 

In  carrying  out  the  computer  calculation,  some  problems  were  encountered 
which  are  worthy  of  note;  in  order  to  establish  a  series  cutoff  criterion, 
it  was  specified  that  when  a  term  contributed  less  than  0. IX  to  the 
preceding  summation,  the  series  was  to  be  terminated.  However,  upon 
attempting  to  verify  this  step  it  was  found  that  it  is  necessary  to  require 
that  not  one,  but  two  successive  terms  be  less  than  0. IX  before  termination 
is  accomplished;  this  requirement  is  due  to  the  oscillatory  nature  of  the 
series  terms. 

An  additional  feature  of  the  f  (k)  computation  which  should  also  be 
noted  is  that  for  large  values  of  k,  f  (k)  becomes  very  small  and  if 
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10? 


accuracy  is  to  be  maintained. ,  care  must  be  taken  to  assure  that  sufficient 
terms  are  included  in  the  series. 

2.  Evaluation  of  f~* { r). 

From  Chapter  V  we  have. 


(SO 

r 


d  k  k 


o 


[ 


1 


mi _ 

- /f(k) 


(5.23) 


In  performing  this  evaluation,  tabulated  values  of  f  (k)  were 
employed.  The  initial  attempt  at  this  computation  involved  the  generation 
of  sufficient  f  (k)  vs  k  points  to  allow  a  smooth  curve  to  be  drawn. 

From  this  curve,  additional  values  of  f  (k)  were  determined  by  means  of 
interpolation;  these  values  along  with  those  actually  generated  by  the 
computer  were  used  in  the  tabulation.  Subsequent  investigation  disclosed 
that  /””* (r)  was  extremely  sensitive  to  the  accuracy  of  the  tabulated  f  (k) 
values.  A  detailed  study  was  th-.n  conducted  to  ascertain  the  number  of 
generated  points  necessary  to  insure  the  desired  accuracy. 

An  additional  Independent  study  was  made  to  evaluate  the  affect  of 
various  &  k  increments  on  the  /""*( r)  computation.  Finally,  the  expression 
for  f  (r)  was  carried  out  to  the  point  where  additional  Integration  incre¬ 
ments  contributed  less  than  1.01. 

It  was  stated  in  the  text  that  the  numerical  work  of  Bowers  was 
considered  incorrect.  In  his  paper he  specified  that  tabulated  values 
of  ker2(i 1)  were  used  in  his  calculations;  the  tables  he  referenced  were: 
Tables  of  Integrals,  by  H.  B.  Dwight,  The  Macmillan  Co.,  New  York,  1949. 
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These  tables  only  provide  ten  points  on  an  f  (k)  curve.  It  has  been  clearly 
established,  during  the  course  of  the  present  work,  that  the  Bowers'  treat¬ 


ment  is  incorrect. 
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ERRATA 


Page  Location 

28  Eq.  (3.  16a) 


Read  exp(- 


Item 

r12  r12 

)  instead  of  exp  (-■jjj'  ) 


z' 


45  Line  1  The  sentence  beginning  at  the  end  of  this  line 

should  read:  Each  term  in  this  expansion  can 
be  represented  graphically  by  a  diagram 
which  consists  of  a  system  of  points  repre¬ 
senting  particles,  and  a  number  of  connect¬ 
ing  lines  each  of  which  corresponds  to  a 
factor  f(r..);  the  total  number  of  connecting 
lines  in  t$e  graph  equals  the  number  of  fac¬ 
tors,  f(r..),  in  the  particular  term  under 
conside  ration. 


52 

Line  18 

3  3  k 

(pr^)  should  read  (pr^  ) 

57 

Line  1 

3 

4tt Gp  =  1  should  read  4itGpX  =  1 

70 

Eq.  (5.20) 

[In  (ek)  -  P(n  +  1)  -  T(n  -  1)]  should  read 
[In  (6k)  -  Y(n  +  1)  -  l(n  -  1)] 

88 

Lines  1 5  and 

17  Interchange  the  words  overestimates  and 
underestimates 

89 

Line  8 

The  sentence  beginning  on  this  line  should 
read:  Therefore,  it  would  be  expected  that 
the  corrected  microfield  distribution  would 
be  shifted  little  if  any,  but  that  the  peak 
would  be  raised  and  the  distribution  nar¬ 
rowed  slightly. 

101 

Eq.  (C.  6) 

Insert  the  factor  — -  in  front  of  the  first 

Vu  + 1) 

integral  inside  brackets;  change  the  upper 
limit  of  integration  in  this  integral  from 
r  to  R. 


